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I. INTRODUCTION 



An important application of lattice gauge theory is to calculate hadronic matrix elements rele- 
vant to experiments in flavor physics. With recent advances in lattice calculations with n/ = 2 + 1 
flavors of dynamical quarks [1-4], we now have an exciting prospect of genuine QCD calcula- 
tions. To match the experimental uncertainty, available now or in the short term, it is essential to 
control all other sources of theoretical uncertainty as well as possible. An attractive target is to 
reduce the uncertainty, from any given source, to 1-2%. This target will be hard to hit if one relies 
on increases in computer power alone: methodological improvements are needed too. 

Many of the important processes are electroweak transitions of heavy charmed or 6-flavored 
quarks. A particular challenge stems from heavy-quark discretization effects, because rriQa ^ 1. 
The key to meeting the challenge is to observe that heavy quarks are non-relativistic in the rest 
frame of the containing hadron [5, 6]. The scale of the heavy-quark mass, rriQ, can (and should) be 
separated from the soft scales inside the hadron and treated with an effective field theory instead 
of computer simulation. Even so, at available lattice spacings [1], many calculations of D-meson 
(£>-meson) properties suffer from a discretization error of around 7% (5%) [2, 3]. Thus, it makes 
sense to develop a more accurate discretization. 

In this paper we extend the accuracy of the "Fermilab" method for heavy quarks [7] to include 
in the lattice action all interactions of dimension six. We also include certain interactions of 
dimension seven. Because heavy quarks are non-relativistic, they are commensurate with related 
dimension-6 terms, in the power counting of heavy-quark effective theory (HQET) for heavy-light 
hadrons [5] or non-relativistic QCD (NRQCD) for quarkonium [6]. 

The Fermilab method starts with Wilson fermions [8] and the clover action [9]. With these 
actions lattice spacing effects are bounded for large mqa, thanks to heavy-quark symmetry. They 
can be reduced systematically by allowing an asymmetry between spatial and temporal interac- 
tions. Asymmetry in the lattice action compensates for the non-relativistic kinematics, enabling 
a relativistic description through the Symanzik effective field theory [10]. Alternatively, one may 
interpret Wilson fermions non-relativistically from the outset [7], and set up the improvement 
program matching lattice gauge theory and continuum QCD to each other through HQET and 
NRQCD [11, 12]. The Symanzik description makes it possible to design a lattice action that be- 
haves smoothly as rriqa — > 0, converging to the universal continuum limit. The HQET description, 
on the other hand, makes semiquantitative estimates of discretization errors more transparent. 

The new action introduced below has nineteen bilinear interactions beyond those of the asym- 
metric version of the clover action, as well as many four-quark interactions. Several of these 
couplings are redundant, and many more vanish when matching to continuum QCD at the tree 
level. We study semiquantitatively how many of the new operators are needed to achieve 1-2% 
accuracy. We find, in the end, that only six new interactions are essential for such accuracy. The 
action is designed with some flexibility, so that one may choose the computationally least costly 
version of the action. 

This paper is organized as follows. Section II considers the description of lattice gauge theory 
via continuum effective field theories. Then, in some detail, we identify a full set of operators de- 
scribing heavy-quark discretization effects. We then determine how many of these are redundant, 
and which redundant directions should be used to preserve the good high-mass behavior. We have 
two goals in this analysis. One is to design the new, more highly improved, action; for this step a 
Symanzik-like description is more helpful, and the resulting action is given in Sec. III. The other 
is to estimate the discretization errors of the new action; here the HQET and NRQCD descrip- 
tions are more useful. To make error estimates, and to use the new action in numerical work, we 
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need matching calculations; they are in Sec. IV. Our error estimates are in Sec. V. Section VI con- 
cludes. Some of the material is technical and appears in appendices: Feynman rules needed for the 
matching calculation are in Appendix A; some details of the Compton scattering amplitude used 
for matching are in Appendix B; a discussion of improvement of the gauge action on anisotropic 
lattices (which one needs only if the heavy quarks are not quenched) is in Appendix C. Some of 
these results have been reported earlier [13]. 



H. EFFECTIVE FIELD THEORY 

In this section we discuss how to understand and control discretization effects using effective 
field theories. We start with a brief overview, focusing on issues that arise for heavy quarks, those 
with mass txiq 3> A. For more details, the reader may consult earlier work [7, 11, 12, 14, 15] 
or a pedagogical review [16]. Here we catalog all interactions of dimension 6 and also certain 
interactions of dimension 7 that, for heavy quarks, are of comparable size when rriQa ^ 1. 



A. Overview 

Cutoff effects in lattice field theories are most elegantly studied with continuum effective field 
theories. The idea originated with Symanzik [10] and was extended to gluons and light quarks by 
Weisz and collaborators [9, 17-19]. One develops a relationship 

Aat = £sym, (2.1) 

where = means that the two Lagrangians generate the same on-shell spectrum and matrix elements. 
The lattice itself regulates the ultraviolet behavior of the underlying (lattice) theory £i at . On the 
other hand, a continuum scheme, which does not need to be specified in detail, regulates (and 
renormalizes) the ultraviolet behavior of the effective theory £s ym - 

In lattice QCD (with Wilson fermions), the local effective Lagrangian (LE£) is 

£ Sym = tr[i^FH - Qf(P+ m f) a f + a^^Kitf, ma; Cj ; ^a)C h (2.2) 

where g 2 and m/ are the gauge coupling and quark mass (of flavor /), renormalized at scale 
fj, < a^ 1 . The (continuum) QCD Lagrangian appears as the first two terms. The sum consists of 
higher dimension operators multiplied by short-distance coefficients K { . These terms describe 
cutoff effects. The short-distance coefficients depend on the renormalization point and on the 
couplings, including couplings Cj of improvement terms in £ lat . Equation (2.2) is fairly well- 
established to all orders in perturbation theory [20, 21] and believed to hold non-perturbatively 
as well. If a is small enough, the terms may be treated as operator insertions, leading to a 
description of lattice gauge theory as "QCD + small corrections". 

In heavy-quark physics uiq 3> A, where A is the QCD scale, so one is led to consider what 
happens when mqa 1. The short-distance coefficients depend explicitly on the mass. Time 
derivatives of heavy-quark or heavy-antiquark fields in the £j also generate mass dependence of 
observables. With field redefinitions — or, equivalently, with the equations of motion — these time 
derivatives can be eliminated. Focusing on a single heavy flavor Q, the result of these manipula- 
tions is [7, 14, 15] 

£ Sym = Q \1aD 4 + mi + -D^JQ + Y, a^^Mg 2 , m 2 a; ^a)£, (2.3) 
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where the ellipsis denotes the unaltered LE£ for gluons and light quarks. By construction the Li 
do not have any time derivatives acting on quarks or antiquarks. 

The advantage of Eq. (2.3) is that all dependence on the heavy-quark mass is in the short- 
distance coefficients mi, \/mi/m 2 , and Ki(m 2 a). Matrix elements of the Li generate soft scales. 
The heavy-quark symmetry of Wilson quarks (with either the Wilson [8] or Sheikholeslami- 
Wohlert [9] actions) guarantees that the coefficients Ki(m 2 a) are bounded for all m 2 a. This fea- 
ture can be preserved by improving the lattice Lagrangian with discretizations of the Li, thereby 
avoiding higher time derivatives [7, 11]. For such improved actions, Eq. (2.3) neatly isolates the 
potentially most serious problem of heavy quarks into the deviation of the coefficient ^/mi/m 2 
from 1. 

Fortunately, the problem can be circumvented in two simple ways. One is a Wilson-like action 
with two hopping parameters [7], tuned so that mi = m 2 . Then Eq. (2.3) once again takes the 
form "QCD + small corrections". The new lattice action introduced in Sec. Ill has two hopping 
parameters for this reason. 

Another solution is to interpret Wilson fermions in a non-relativistic framework. One can 
replace the Symanzik description with one using a non-relativistic effective field theory for the 
quarks (and antiquarks) [1 1]. For the leading Q-Q term in Eq. (2.3) 

QL Di + m 1+ ./™T 7 .d)o = SW (D l + mi - D2 + ZB(m : a ^ a) ^- B )h^ + ... 
\ ]/ m 2 J V 2m 2 J 

(2.4) 

where z B is a matching coefficient, and is a heavy-quark field satisfying = +'j i h ( - + \ 
Another set of terms appears for the antiquark, with field h^' satisfying hS' = — 74M - ). The 
non-relativistic effective theory conserves heavy quarks and heavy antiquarks separately. As a 
consequence, the rest mass m 1 has no effect on mass splittings and matrix elements. 1 For lattice 
gauge theory this implies that the bare quark mass (or hopping parameter) should not be adjusted 
via mi. Instead, the bare mass should be adjusted to normalize the kinetic energy D 2 /2m 2 . 

One can develop the non-relativistic effective theory for the lattice artifacts Li by using heavy- 
quark fields instead of Dirac quark fields [11]. Higher-dimension operators in the heavy-quark 
theory receive contributions from the expansions of Eq. (2.4) and of the Li. Coalescing the coeffi- 
cients of like operators obtains a description of lattice gauge theory with heavy quarks 

Aat = • • ■ - U+\D 4 + mi)h^ +J2Cl at {9 2 ,rn 2 ;m 2 a,c j - f i/m 2 )O i , (2.5) 

i 

where the operators Oi on the right-hand side are those of a (continuum) heavy-quark effective 
theory, of dimension 5 and higher, built out of heavy-quark fields h^', gluons, and light quarks. 
(The leading ellipsis denotes term for the gluons and light quarks only.) The C, are short-distance 
coefficients, which depend on g 2 , the heavy-quark mass, the ratio of short distances m 2 a, and also 
all couplings cj in the lattice action. The logic and structure is the same as the non-relativistic 
description of QCD, 

£qcd = • • • - h {+) (D 4 + m Q )h^ + ]T CrV, rn Q ; /i/m Q )0, (2.6) 

i 

Thus, improvement of lattice gauge theory is attained by adjusting couplings Cj until C] at (c,,) — Cf ont 
vanishes (identically, or perhaps to some accuracy) for the first several Oi. 



A simple proof can be found in Ref. [11]. 
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TABLE I: Bilinear interactions that could appear in the Symanzik LEC through dimension 6. 
Dim With axis-interchange symmetry Without axis-interchange symmetry HQET A s NRQCD v t 

3 qq QQ 



4 qpq 


0(741)4 + mi)Q 




1 


v 2 




Ql-DQ 




A 


v 2 


5 q~D 2 q ei 


QDjQ 


£1 








QD 2 Q 


<h 


A 


v 2 


— TiPUmjFinjQ 


QiYt ■ BQ 




A 


v 4 




Qa ■ EQ 




A 2 


v 4 


6 qjuD^q 
" if* fiJ- 


QliDfQ 




A 3 


v 4 


q{B,D 2 \q e 2 


QiaDIq 


£2 








Ql7 4 -D 4 , -D 2 |Q 


6 2 








g{D|, 7 .u}g 


■d 2 








Q{ 7 -D,D 2 }Q 




A 3 


v 4 


-^q{Ip,a^F^}q e F 


Q{ 7 - D,a £}Q 


£ F 


A 2 


v 4 




Q{ 74 D 4 ,^-S}Q 


Sb 








Q{ 7 -D,tS-B}Q 




A 3 


V 6 




Q[£»4, 7 --E]Q 




A 3 


V 6 


q[D^,F^ u q 


074(0 E-E D)Q 




A 2 


V 4 




Q-f-(DxB + Bx D)Q 




A 3 


V 6 



It does not matter whether one carries out the improvement program by adjusting Ki(cj) = 
or C] at (cj) = Ci° nt [12]. The results for the c,- are the same, provided one identifies txiq with m 2 . 
The matching assumes that pa <C 1, but at the same time m 2 a 1. One is thus led to 
non-relativistic kinematics (p/m 2 <C 1) in the matching calculation, where both descriptions — 
Eqs. (2.3) and (2.5) — are valid. Kinematics are encoded into the operators £; or Oi and are not 
transferred to the short-distance coefficients. Hence, kinematics cannot influence matching condi- 
tions on the Cj. In particular, when indeed m 2 a <C 1 (which may be impractical, but is conceivable 
theoretically) relativistic kinematics (p ~ m 2 ) are possible, and it follows from the Symanzik 
effective field theory that the solution of K^Cj) = yields the same cj for both relativistic and 
non-relativistic kinematics. 

B. Quark bilinears in the LE£ 

In the rest of this section we construct the LE£ appropriate to heavy quarks. The two main 
steps are first to list all of the Cj that can appear, and second to decide which should be considered 
redundant. In part it is a generalization of the dimension-6 analysis of Ref. [9] to the case without 
axis -interchange symmetry. At dimension 6 there are quark bilinears, four-quark interactions, and 
interactions that contain only the gauge field. We shall start with the bilinears and turn to the others 
further below. In each case, we first consider complete lists of operators, and then consider which 
can be chosen to be redundant. 

Table I contains a list of all quark bilinears through dimension 6 that can appear in the effective 
Lagrangian. The second column contains interactions that respect axis-interchange symmetry; the 
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fourth column contains the extension to the case without axis -interchange symmetry. The meaning 
of the other columns is explained below. Covariant derivatives act on all fields to the right, 



D^FQ = (d^F + [Ap, F])Q + F D^Q. (2.7) 

This notation is convenient for the interactions with commutators and anti-commutators. To arrive 
at the lists we exploit identities such as 

f = D 2 - i<vF^, (2.8) 
274D47 • £>7 4 £ 4 = {74^4, ot-E}- {Dj, 7 ■ D}, (2.9) 
2 7 • D74D47 D = {7 D, a -E} - {74D4, (7 • Df}. (2.10) 



Some interactions are omitted, because the underlying lattice gauge theory is invariant under cubic 
rotations, spatial inversion, time reflection, and charge conjugation. 2 

The fourth column is arranged so that its entries are part of the corresponding interactions in the 
second column. It is easy to show that the list is complete, by writing out all independent ways to 
have three covariant derivatives, expressing the E and B fields as anti-commutators of covariant 
derivatives. One finds 1 1 possibilities, and then one can use identities to manipulate this list to that 
given in the fourth column of Table I. 

The LE£ contains several redundant directions. The equation of motion of the leading LE£ 
plays a key role in specifying which operator insertions may be considered redundant. Let us 
assume, for the moment, that m\ = m 2 , so that the equation of motion in the Symanzik LE£ is the 
Dirac equation. Below we shall use the non-relativistic effective field theory to address the case 

m l 7^ m 2- 

The quark fields are integration variables in a functional integral, so an equally valid description 
is obtained by changing variables 

Q >-> e J Q, (2.11) 

Q >-> Qe J , (2.12) 

where 

J = ae 1 (p+ m) + a5i7 • D + a 2 e 2 {Ip + m) 2 - a 2 \e F io ^F^ + a 2 5 2 (-y ■ D) 2 

+ a 2 S B iYl ■ B + a 2 tf 2 [7 4J D 4 ,7 • D] (2.13) 

and similarly for J with separate parameters e i} Si, and If the 5 parameters (and $2, ^2) vanish, 
then J and J preserve invariance under interchange of all four axes. 

One can propagate the change of variables to the LE£, and trace which coefficients of dimen- 
sions 5 and 6 are shifted by amounts proportional to the parameters in J and J. To avoid generating 
terms that violate charge conjugation one chooses E{ = +£i, 5i = +Su ^2 = —$2- We then see that 
there are two redundant directions at dimension 5, and five at dimension 6. That means that two 
couplings in the dimension-5 lattice action may be set by convenience, and five in the dimension-6 
lattice action. The third and fifth columns show the correspondence between parameters in the 
change of variables and the interactions that we choose to be redundant. As expected from general 



2 Reference [9] included the dimension-6 interaction q[p, D 2 ]q. Reference [7] included the dimension-5 interaction 
Q[74£ ) 4, 7 • D]Q- Both are odd under charge conjugation and, thus, may be omitted. 
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arguments [7, 14, 15], all interactions in which 74-D4 acts on Q or (after integration by parts) Q 
are redundant. 

There is quite a bit of freedom here. One could choose Ep to eliminate Q[D^~{ ■ E]Q = 
Q{^aD Ai ol ■ E}Q instead of Q{j ■ D,az ■ E}Q. But the former is suppressed, relative to the 
latter, in heavy-quark systems. Moreover, in HQET and NRQCD one has 

Qa-EQ = h {+) { 1 -D,a-E}h {+) /2m 2 + ---, (2.14) 
Q{7 D,a E}Q = h {+) {>y ■ D, a ■ E}h {+) + ■ • • , (2.15) 

which mean that Qct ■ EQ and Q{j ■£>,«■ E}Q generate nearly the same effects in heavy-quark 
systems. Thus, we prefer to take Q{j ■ D, a. ■ E}Q to be redundant. 

To understand the general pattern of redundant interactions, let us introduce some notation. 
Let B (£) be a combination of gauge fields, derivatives, and Dirac matrices that commutes (anti- 
commutes) with 74. An example of B (£) is zS ■ B (a ■ E). Also, let us write B± (and £±) when 
QB±Q (or Q£±Q) has charge conjugation ±1. Because we wish to eliminate time derivatives 
of quark and antiquark fields, we would like Q^^D^, B + }Q and Q[7 4 -D 4 , £-}Q to be redundant. 
That is always possible: simply add to J in Eq. (2.13) terms of the form 5b + B + and As 
a consequence, neither Q{j ■ D, B + }Q nor Q[~f ■ D, £-}Q is redundant. On the other hand, in 
Ql'jiDi, B-]Q and Q{74-D 4 , £+}Q the time derivative acts only on gauge fields. Thus, by adding 
to J terms of the form tDqB^ and §e + £+ it is possible to choose Q [7- D, BJ\Q and QIj-D, £ + }Q 
to be redundant. Instead of Q[j ■ D, BJ\Q or Q{j ■ D, £ + }Q it may be convenient to choose an 
operator related through an identity. 

C. Power counting 

The small corrections of an effective field theory are small, because the product of the short- 
distance coefficients and the operators yield a ratio of a short-distance scale to a long-distance 
scale. For light quarks in the Symanzik effective field theory, the essential ratio is a/A -1 = Aa, 
and dimensional analysis reveals the power of Aa to which any contribution is suppressed. In 
particular, B- and £-type interactions of the same dimension are equally important. 

For heavy quarks the physics is different, because m^ 1 is a short distance. The ratio a /nig 1 = 
mQd should not be taken commensurate with Aa [7]. Instead, interactions should be classified in 
a way that brings out the physics. It is natural to turn to HQET and NRQCD. Let us start with 
heavy-light hadrons and HQET. £-type interactions of given dimension are A /nig times smaller 
than £>-type interactions of the same dimension. Because A/uiq 1 and Aa 1, it makes sense 
to count powers of A, where A is either of the small parameters [11, 12, 15] 

A ~ aA, A/rriQ. (2.16) 

This power counting pertains whether rriq < a, vtlq ~ a, or tuq > a. Writing the corrections in 
the Symanzik fashion (with Dirac quark fields Q and Q), each Li is suppressed by A s , with 

s = dim£-4 + n r . (2.17) 

Here n r = or 1 for interactions of the form QB + Q or Q£ + Q, respectively. The sixth column 
of Table I (labelled HQET) shows the suppression of each interaction, relative to the (leading) 
contribution from the light degrees of freedom. In the following we call the power counting for 
heavy-light hadrons, based on Eq. (2.17), "HQET power counting." 
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TABLE II: Dimension-(7,0) bilinear interactions that are commensurate, for heavy quarks, with those of 
order A 3 (in HQET) or v 4 , v 6 (in NRQCD). 



Without axis-interchange symmetry 


HQET A s 


NRQCD 


QDfQ 




A 3 


v 4 


J2i^j QiXiDjBiDjQ 




A 3 


V 6 


Ei^Q{D],i^iBi}Q 




A 3 


V 6 


Q{D 2 ) 2 Q 




A 3 


V 1 


Q{D 2 ,iY,-B}Q 




A 3 


V 6 


Q7 • DiY, • B7 • DQ 


5[{ 7 . 


A 3 


I > 6 


QDiiS ■ BDiQ 




A 3 


I > 6 


QD ■ (B x D)Q 


5 [7 • (D x B + B x D)] 


A 3 


v 6 


Q(i£ • B) 2 Q 


«5[{7"A-D 2 }] 


A 3 


v 8 


QB-BQ 




A 3 


v s 


Q{cx-E) 2 Q 


ap>4,7--E]] 


A 3 


V 6 


QE-EQ 




A 3 


v 6 



Now let us recall how to classify interactions in quarkonium according to the power of the rel- 
ative internal velocity, v. Because color source and sink are both non-relativistic, chromoelectric 
fields carry a power of v 3 , and chromomagnetic fields a power of v 4 [22]. £-type interactions are 
suppressed by a power of vl m Q — v > analogously to their suppression in heavy-light hadrons. 
Thus, bilinears are suppressed by v*, where now 

t = dim£ — 3 + n E + 2n B + rip, (2.18) 

and tie (ns) is the number of chromoelectric (chromomagnetic) fields. The seventh column of 
Table I (labelled NRQCD) shows the suppression of each interaction. In the following we call the 
power counting for quarkonium, based on Eq. (2.18), "NRQCD power counting." 

Glancing down the sixth and seventh column of Table I, one sees several terms of order A 3 
and v 6 , from Eqs. (2.17) and (2.18) one realizes that some dimension-7 interactions are of the 
same order. They are listed in Table II. There are two interactions with four derivatives, six with 
the chromomagnetic field and two derivatives, and four with two E or two B fields. A third 
combination of four derivatives is omitted, using the identity DiD 2 Di = (D 2 ) 2 + D ■ (B x D) — 
B 2 . Other dimension-7 operators carry power A 4 in HQET power counting, or t; 8 (or higher) in 
NRQCD power counting. Five combinations are redundant (as shown), and we shall see below 
how they and the others arise in matching calculations. 

The (d, nr) = (7, 1) operator Q{D 2 , a ■ E}Q and several (d, n r ) = (8, 0) operators, all with 
tie — 1 and «o + nr = 3, have NRQCD power-counting v 6 . Reference [22] includes spin- 
dependent ones, to obtain the next-to-leading corrections to spin-dependent mass splittings. We 
have not included these operators in our analysis, but a straightforward extension of the matching 
calculation in Sec. IV B 1 would suffice to determine their couplings. 

Although this description of cutoff effects is somewhat cumbersome, it provides a valuable 
foundation for our new action, given in Sec. III. To obtain the new action, we simply discretize 
the interactions in Tables I and II, except those with higher time derivatives. The discretization of 
Qj ■ DQ is needed to obtain a lattice action that behaves smoothly as rriQa — » [7], reproducing 
the universal continuum limit of QCD. Similarly, discretizations of the £-type interactions, such 
as Qct ■ EQ and Q{ 7 ■ D, D 2 }Q, are needed to retain that feature here. 
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TABLE III: Bilinear interactions that could appear in the heavy-quark LEC through dimension 7. 


Dim Without axis-interchange symmetry 


HQET A s 


NRQCD v* 


3 hWh& 

4 /i( ± ) 74J D 4 /i (±) 


5 h&Dlh& £i 








A 


v 2 




A 


V 4 


6 ^( ±) 7 4 Z?|/i (±) e 2 






/j (±) {74-D4,-D 2 }/i (±) 5 2 






/j(±){ 7 D,a ■ E}h^ 


A 2 


V 4 


^ (±) {74D 4 ,iS • Sj/tW <5b 






P T4 (D -E-E- D)h^ 


A 2 


v 4 


1 h&Dfh& 


A 3 


v 4 




A 3 


v % 


^■h^i^DjBiDjh^ 


A 3 


V G 


m{D 2 ) 2 h^ 


A 3 


V 4 


h^{D 2 ,iS- B}h^ 


A 3 


,fi 


h^j ■ DiT, • B 7 ■ DfcW 


A 3 




U^DiiS ■ BDihW 


A 3 




• (S x D)h^ 


A 3 


V 6 


U^(iV-B) 2 h& 


A 3 


i' 8 


h&B ■ Bh& 


A 3 


v 8 




A 3 


V 6 


hP>E ■ Eh& 


A 3 


V 6 


D. Heavy-quark description 

For understanding the size of heavy-quark discretization effects, it is simpler to switch to a 



non-relativistic description. (When m 1 ^ m 2 , it is also necessary to see the connection to QCD.) 
The list of interactions is much shorter, because the constraint "f^h^ = ±h^ removes the S- 
type interactions. It is given in Table III, including the dimension-7 interactions related to those in 
Table II. Also, fewer changes of the field variables are possible: 

h {±) i — ► e J h, (2.19) 

h(±) i — > he J , (2.20) 

where now 

J = ae 1 (-f 4 D 4 + mi) + a 2 £ 2 (74£>4 + »i) 2 + a 2 5 2 D 2 + a 2 5 B z'E • B, (2.21) 

and similarly for J. To avoid C-odd interactions, one should choose equal parameters in J and J. 
Thus, there are four redundant directions of interest — all with time derivatives of the (anti-)quark 
field. In the end, just as many non-redundant interactions remain as in the Symanzik description. 
The heavy-quark description provides a good way to estimate the size of remaining discretization 
effects, as in Sec. V. 



9 



E. Gauge-field and four-quark interactions in the LE£ 

We now turn to interactions in the gauge sector of the LE£, and also to four-quark interactions. 
The two are connected when one considers on-shell improvement, because in quark-quark scatter- 
ing short-distance gluon exchange generates the same behavior as four-quark contact interactions. 
Here we give a cursory sketch of the gauge action. Then we consider the four-quark interactions, 
including details mostly for completeness. In practice (see Sec. V), we find the four-quark correc- 
tions to be smaller than those of the bilinear interactions analyzed in the preceding subsection. 

The gauge sector of the LE£ is the same as for anisotropic lattices, where one adjusts the action 
so that the temporal lattice spacing a t differs from the spatial lattice spacing a s . The short-distance 
coefficients are different; here asymmetry between spatial and temporal gauge couplings arise only 
from heavy-quark loops. Improved anisotropic actions have been discussed in the literature [23], 
but full details remain unpublished [24]. We present the details in Appendix C. 

We are most concerned here with effects that survive on shell, so we study here the possible 
changes of variables for the gauge field. With axis -interchange symmetry one has [9, 19] 

A^A, + a 2 e A [D\ F^} + a 2 g 2 ]T e Jf t a (q f ^t a q f ), (2.22) 

/ 

with a color-adjoint vector-current term for each flavor / of quark (heavy or light). The appearance 
of g 2 multiplying the currents is a convenient normalization convention. When one now considers 
giving up axis -interchange symmetry, one has 

A A h-f A A + a 2 e A (D ■ E - E ■ D) + a 2 g 2 e Jf t a (q f ^t a q f ), (2.23) 

/ 

A i-> A - a 2 (e A + 5 E )[D A , E] + a 2 (e A + 5 A )(D x B + B x D) 

+ a 2 g 2 J2(tJf + S Jf )t a (q n t a q f ), (2.24) 
/ 

which reduce to Eq. (2.22) when the 5s vanish. 

For a moment, let us set ejj = Sjf = in Eqs. (2.23) and (2.24), and focus on the gauge 
fields alone. As discussed in Appendix C, there are eight independent gauge-field interactions that 
arise at dimension six. There are three independent ways — parametrized by e A , 5 A , and 5e — to 
transform the gauge field, yielding three redundant directions. Similarly, there are eight distinct 
classes of six-link loops, shown in Fig. 1, that can be used in an improved lattice gauge action. In 
Appendix C, we show that three of them — all three classes of "bent rectangles" in the bottom row 
of Fig. 1 — may be omitted from an on-shell improved gauge action. 

The transformations involving the currents 9/7 M t a ?/ are more interesting. They shift the LE£ 
[cf. Eq. (2.2)] by 

^ £sym — fl2 

ej f q f y 4 (D E-E D)q f + a 2 ^{ejf + S Jf )q f [D 4 , 7 • E}q f 
f f 

- a 2 Y^( £ Jf + 8jf)qn -(DxB + Bx D)q f 

f 

' <»VE e J/(?/^9/)(9 ff 7/?j) -aVXJ^/(^* a 9/)(^%), (2-25) 
fg fg,j 
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FIG. 1: Six-link loops available for improving the gauge action on anisotropic lattices: rectangles (top row); 
parallelograms (middle); bent rectangles (bottom). Nomenclature from Ref. [19]. 

where the derivatives act only on the gauge fields. The size of these shifts — of order g 2 for four- 
quark operators and of order g° for bilinears — is commensurate with the respective terms that al- 
ready appear in £s ym - Thus, the 2nj parameters ejf and 5jf could be used to eliminate bilinears or 
four-quark operators. For simulations it is preferable to remove the latter, namely qfJi^QfQfJi^qf 
and qf~lt a qf ■ qf~it a qf. 

We now list the dimension-six four-quark interactions in the LE£. For a single flavor, the 
complete list is in Table IV, which also indicates that the current-current interactions are redundant. 
Interactions with the color structure (qTq) 2 may be omitted, because they can be related to those 
listed through Fierz rearrangement of the fields. 

When considering several flavors of quark, we must keep track of flavor indices as well as color 
and Dirac indices. The Fierz problem becomes more intricate, and we shall find that color- singlet 
and color-octet structures should be maintained. Let us start with Fierz rearrangement of the Dirac 
indices. The four-quark terms in the LE£ take the form 



KxQf a Tx<lg/3Qhj^X<liS = - ^ KxFxYqfa^YqiSQhy^YqgP, 



(2.26) 



X.Y 



where K x denotes short-distance coefficients, the Greek (Latin) indices label color (flavor), F is 
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TABLE IV: Four-quark interactions that could appear in the LEC (for a single flavor). 



Dim With axis interchange Without axis interchange 

~6 (qt a q) 2 (Qt a Q) 2 

(qi5t a q) 2 {Ql^Q) 2 
(q^t a q) 2 ej (QlAt a Q) 2 ej 

(Qlit a Q) 2 Sj 
(qi^t a q) 2 (Q lAl5 t a Q) 2 

(Q^t a Q) 2 

{qia^t a q? (Qi^it a Q) 2 



{Qait a Q) 



2 



TABLE V: Four-quark interactions that remain when Fierz rearrangement is taken into account. A sum over 
Dirac matrices T x in each of the sets {1}, {74}, {7}, {i^}, {ex}, {775}, {7475}, {75} is assumed. (With 
axis-interchange symmetry, the sets would be {!}, {7^}, {icr^ v }, {7^75}, {75}-) 



Quarks 


Color octet 


Color singlet 


Heavy-heavy 


QT x t a QQT x t a Q 




Heavy-heavy 


QiT x t a Q l Q 2 T x t a Q 2 


QiT x QiQiTxQ2 


Heavy-light 


QT x t a Q^ f q f T x t a q f 


QTxQT,fqf T xqf 


Light-light 


T. f q^xt a q f T. q q^xt a q g 





the Fierz rearrangement matrix (with F 2 = 1), and the minus sign comes from anti-commutation 
of the fermion fields. Equation (2.26) leaves the flavor and color indices uncontracted, but to get 
terms in the LE£, the color indices must be contracted (one way or another), and the flavor labels 
must yield a flavor-neutral interaction. Without loss, we can choose the side of Eq. (2.26) such 
that the Dirac matrices contract quark fields of the same flavor. Then one can use Fierz identities 
for SU(N) generators (t at = -t a ) 

Nt a a fy = -t a aS t a y0 - (N 2 - l)5 aS 5^/2N, (2.27) 
<W<5 7<5 = 6 aS 6 7 p/N-2t a aS t^, (2.28) 

so that the color indices are contracted across the same fields as the Dirac and flavor indices. 

After using Fierz rearrangement to bring quarks of the same flavor next to each other, one 
is left with the interactions in Table V. To be concrete, we consider n; flavors of light quarks 
(with m q < A) and two flavors of heavy quarks (charm and bottom). We neglect the dependence 
of the coefficients on the light quark masses, because four-quark interactions are already small 
corrections (of dimension six). In that case, the four-quark interactions can be arranged so that 
only the SU(n;) flavor singlets J2f Qf^xt a Qf an d qfT x qj appear. 

The parameters sjf and Sjf may be used to eliminate color-octet current-current interactions. 
For each heavy flavor, one finds (Qj^Q) 2 and J^iiQ'li^Q) 2 to be redundant. For light quarks, 
we may neglect the differences in the mass, so they have common parameters, and the flavor- 
singlet combination (J2f Qfl^Qf) 2 ls redundant. For the light flavors, our list of operators is a 
Fierz rearrangement of the list in Ref. [9]. 

The leading HQET power counting for heavy-light four-quark operators follows from dimen- 
sional analysis and Eq. (2.17): A 2+nr , just as if the light-quark part were replaced by three deriva- 
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tives. Heavy-heavy four-quark operators will be suppressed, once matrix elements are taken, by a 
heavy-quark loop, leading to g 2 \ 4+nv . 

In quarkonium, the size of heavy-light four-quark operators follows similarly from Eq. (2.18): 
v 3+nr . The valence heavy-heavy operators are more interesting. They must contain two contribu- 
tions, one to improve t-channel gluon exchange, and another to improve s-channel annihilation. 
The former have NRQCD power counting g 2 v 3+nr ~ v A+n r (since g 2 ~ v [22]). The latter are v 2 
times smaller, because the s-channel gluon is far off shell, but the Dirac-matrix suppression is now 
f 1_nr , leading to g 2 v 6 ~ nr ~ v 7 ~ nr in all. In practice, the s-channel contributions are suppressed 
further, when treated as an insertion in a color-singlet quarkonium state. At the tree level, the 
only color structure that can arise is the color-octet. Its matrix elements vanish in the QQ-co\or- 
singlet Fock state of quarkonium, leaving the t; 3 -suppressed QQA color octet [25]. Color-singlet 
four-quark operators arise at one loop, with an additional factor of g 2 ~ v. 



HI. NEW LATTICE ACTION 

In this section we introduce a new, improved lattice action for heavy quarks, designed to yield 
smaller discretization errors than the action in Ref. [7]. Our design is based on several lessons 
from the preceding section and Refs. [7, 11, 12]. First, it is important to preserve the natural 
heavy-quark symmetry of Wilson fermions, so that the coefficients Ki stay bounded for all vriQa. 
(This feature is spoiled in the standard improvement program designed for light quarks, which 
introduces several new terms that grow with rriQ.) Second, the new lattice action is flexible enough 
to match cleanly onto both the Symanzik description and the non-relativistic description. 

Let us write the action as follows 

oo 1 

S = S D 2 F 2 + So + 2^ ^ S(d,n r ) + ^9999' (3.1) 
d=5 n r =0 

where S D 2 F 2 is the improved gauge action [Eq. (C7)], S is the basic Fermilab action, the S(d,n T ) 
consist of the bilinear terms added to improve the quark sector, and Sg q g q denotes four-quark 
interactions. S(d,n r ) consists of (discretizations of) interactions of dimension d, with nr as in 
the discussion of power counting, Eqs. (2.16)— (2.18). Including the interactions in S(d,i) couples 
"upper" and "lower" components, but allows a smooth limit a — > 0. 3 Our aim is to improve the 
action to include all interactions of dimension six. Then the power counting requires us to include 
Spfi) as well. Finally, S qqqq consists of discretizations of four-quark operators, at dimension six, 
those of Table V. 

The basic Fermilab action [7] is a generalization of the Wilson action [8]: 

S = m a 4 ^ $(x)ij>(x) + a 4 ^ ^( x h^ D ^t^(x) - |a 5 tff(x) A 4lat ^(x) 

XX X 

+ Ca 4 Y,^h ■ Aat^W - ir s Ca 5 ^^(x)A{ a 3 t V(x). (3.2) 

X X 

We denote lattice fermions fields with ij) to distinguish them from the continuum quark fields 
in Sec. II. The dimension-five Wilson terms are included in 5*0 to remove doubler states. The 



3 Lattice NRQCD, which directly discretizes the continuum heavy-quark action, can be thought of as omitting S^i) 
in favor of S {d+lfl) . 
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remaining dimension-five interactions are [7, 9] 

5(5,0) = S B = -ic B Ca 5 ^^(a;)zE ■ Bi at tp(x), (3.3) 

X 

5(5,1) = S E = -lc E (a 5 J2iP(x)a ■ Ete^x), (3.4) 

X 

(3) 

where the notation S E and S E is from Ref. [7], and the discretizations D nt , A w ., A lat , Aat> 
Aat are defined below. 

The new interactions in Eq. (3.1) introduced in this paper are 

5(6,o) = r E a 6 ^ $i x ){l • Aat, ol ■ AatMx) 

x 

+ Z E a 6 ^(^)74 (Aat • Slat - ^at ' Aat) lj}{x\ (3.5) 
X 

5(6,i) = cia 6 ^^(x)^7iAiatA ilat ^(x) + c 2 a 6 ^^(x){7- Aat, A[2}^(x) 

x i x 

+ c 3 a 6 ^^(x){7- Aat,«S- Aat}^(x) 

x 

+ Z 3 a 6 ^2tp(x)~f ■ (Aat X Aat + Aat X Aat) 4>(x) 
x 

+ c EE a G ${x) {74 Aiat, " ■ Aat}^(x), (3-6) 
5(7,o) = c 4 a 7 XI A ^( x ) + c 5 a 7 ^ rj}[x) ^ Alat, A Jla J^(x) 

x i x i j^i 

x i j^i 

+ 4a 7 £^(x)[A«£ • BAliat^W 

+ r 7 a 7 X ^(x)7 • Aat^S • Aat 7 * AatVK" 

X 

+ r' 7 a 7 ^$(x)[D-(B x ^)]iat^(ar) 

x 

+ r B Ba 7 ^ ^(x) (z'S • Aat) 2 + ^BBO 7 X ^(x) Aat ' Aat^( 3 

x x 

- r EE a 7 X i>(x) (a ■ Aat) 2 ^(s) + z EE a 7 ^ ^(ar) Aat ■ Aat^(z)- (3.7) 

X X 

All couplings in Eqs. (3.2)-(3.7) are real; explicit factors of i arc fixed by reflection positivity [26] 
of the continuum action. Some of the improvement terms extend over more than one timeslice, so 
there are small violations of reflection positivity for the lattice action. We expect that the associated 
problems are not severe, as with the improved gauge action [27]. 
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Equations (3.5)-(3.7) contain 19 new couplings. The convention for couplings q, and Zi 
is as follows. In matching calculations we find that couplings z,- L vanish at the tree level, while 
the couplings q do not. Couplings r, are redundant and, for this reason, could be omitted. The 
analysis in Sect. II gives the number of redundant interactions, rather than the specific choices of 
interactions themselves. The possibilities for the dimension-7 redundant directions are as follows. 
One of (c 4 , c 5 , r 5 ) is redundant; we choose r 5 . Furthermore, one of (z 6 , z 7 , r 7 , r BB ), another of 
(z 7 ,r 7 ,r BB ), and another of (z 7 ,r 7 ,r' 7 ,r BB ) are redundant; we choose r 7 , r' 7 , and r BB . But be- 
cause pragmatic considerations could motivate other choices, we keep all of them in our analysis. 
This strategy also provides a good way for the matching calculations to verify the formal analysis 
of the LE£. In future numerical work, we recommend choosing r s , as usual, to solve the doubling 
problem (in practice r s > 1). The others may be chosen to save computer time, which presumably 
means choosing the couplings of computationally demanding interactions to vanish. 

The difference operators and fields with the subscript "lat" are taken to be 

D nat = (T p -T^)/2a (3.8) 

3 

A Plat = (T P + T„ p - 2)/a 2 , A2 = A *iat, (3-9) 

i=l 

Fpa Ut = S s g n P s S na " [TpTvT_ p T_v - T- a TpT_vT_- p \ , (3.10) 

p=±p cr=±a 

where the co variant translation operators T± p translate all fields to the right one site in the ±p 
direction, and multiply by the appropriate link matrix [28]. These discretizations are conventional 
for S + S B + S E . For the new interactions, we have re-used the same ingredients. 
For the interactions with couplings r 5 and z' 7 one can consider 

[DjBiD^ = D jlat B^ t D j}&v (3.11) 

or 

[DjBiDj]^ = ^ [(1 - T-j^jTj - 1) + (Tj - l)B ilat (l - T.j)] . (3.12) 

In tree-level matching calculation, both lead to the same dependence on r 5 and z' 7 . Equation (3.11) 
has the advantage that is re-uses elements that are already defined (in a computer program, say) for 
the dimension-4 and -5 action. Equation (3.12) is more local, however, and may have other advan- 
tages. A FermiQCD [29] computer code of the new action indicates that Eq. (3.11) is faster [30]. 
This code also indicates that it is advantageous to choose the redundant directions so that one may 
set r 5 = r 7 = 0. 

The improved gluon action S D 2 F 2 is defined in Appendix C. The four-quark action Sq q q q con- 
tains the obvious discretization of the (continuum) operators explained in Sec. HE and listed in 
Tables IV and V: simply substitute lattice fermion fields for the continuum fields, and assign 
each a real coupling. When matching to continuum QCD, the couplings in Sq qq - q start at order g 2 , 
making them commensurate with order-g 2 matching effects in SV 6)1 ) + S/ 7i0 ), such as tree-level 
quark-quark scattering. To incorporate the four-quark action in a Monte Carlo simulation, one 
would introduce auxiliary fields to recover a bilinear action. In the next section we show, however, 
that these operators are not necessary for the target accuracy of 1-2%, so this cumbersome set-up 
can be avoided for now. 
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IV. MATCHING CONDITIONS 



In this section we derive improvement conditions on the new couplings at the tree level. We 
calculate on-shell observables for small pa without any assumption on mqa. We look at the 
energy as a function of 3-momentum, which is sensitive to c\, c 2 , c 4 , and z Q . We then look at 
the interaction of a quark with classical background chromoelectric and chromomagnetic fields. 
The former is sensitive to c E , r E , and z E ; the latter to all but c EE , r EE , z EE , r BB , and z BB . To 
ensure that these results are compatible with the improved gauge action, we next compute the 
amplitude for quark-quark scattering. This step also matches the four-quark interactions, which 
are not written out explicitly in Sec. III. Finally, we compute the amplitude for Compton scattering 
to match c EE , r EE , z EE , r BB , and z BB . 



A. Energy 

The energy of a heavy quark on the lattice is defined through the exponential fall-off in time of 
the propagator. For small momentum p the energy can be written 



E = nil + 



PL 

2m 2 



UP 



+ 



5m, 



(4.1) 



where the coefficients mi,m 2 , m 4 and u> 4 depend on the couplings in the action. Appendix A con- 
tains the Feynman rule for the propagator and recalls the general formula for the energy, Eq. (A4). 
By explicit calculation we find 



mid 
1 

m 2 a 
1 

m?a 3 



ln(l + m a), 
2( 2 

m a{2 + m a) 
2C(C + 6 Cl ) 



+ 



1 + m a 
r s ( - 24c 4 



moa(2 + rriQa) 4(1 + rrioa) 
' 8C 4 



(4.2) 
(4.3) 

(4.4) 



[moa(2 + m a)] 3 
| 32(c 2 
m a(2 + m Q a) 1 + mod 



| 4C 4 + 8r s C 3 (l + m a) | 
[m a(2 + m a)} 2 
8z 6 



1 + m a) 2 



(4.5) 



The dimension-6 and -7 couplings (ci, c 4 ) and (c 2 , z G ) modify u> 4 and m 4 a, but not m 4 a or m 2 a. 

To match Eq. (4.1) to the continuum QCD, one requires m 4 = m 2 and u> 4 = 0. From m 4 = m 2 
one obtains the tuning condition 



4C 4 (C 2 



16Cc 2 



C 3 [2C + 4r s (l + m a) - 6r s ( 2 /(l + m a)\ 



[m a(2 + m a)]^ 



2a4 



+ 



3ri( 



m a) 



+ 



m a(2 + m a) 
2(1 + mod) 



8^6 + 



rriQa(2 + m^a) 

^3/^ 3 



2a2 



moa) ; 



mod 



(4.6) 



which (at fixed m a) prescribes a line in the (c 2 , z e ) plane. From W4 = one obtains the tuning 
condition 

-2 , c/ - r„ /• o/i \ m o a (2 + m a) 



= C + Kci + (r s ( - 24c 4 ) 



8(1 + m a) 



(4.7) 
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which (at fixed m a) prescribes a line in the (ci,c 4 ) plane. As m a — > 0, both lines become 
vertical: the coefficients c\ and C2 of dimension-6 operators are fixed, whereas the coefficients of 
C4 and zq dimension-7 operators are undetermined. At this stage it is tempting to choose C4 and z 6 
to be two of the redundant couplings, but below we shall see that there are better choices. 



B. Background Field 



To compute the interaction of a lattice quark with a continuum background field, we have to 
compute vertex diagrams with one gluon attached to the quark line. The Feynman rules are given 
in Eqs. (A23) and (A24). Our Feynman rules introduce a gauge potential via 

U^x) = exp [goA^x + |e M a)] , (4.8) 

where is a unit vector in the fi direction, and take the Fourier transform of the gauge field to be 

M*) = j (0 e * fc ' x A^)- (4-9) 

A background field would, however, lead to parallel transporters 



U^x) = Pexp 



go / A^x + se^ajds 



(4.10) 



Equation (4.8) is a convention. If we use Eq. (4.10) instead, vertices, propagators, and external line 
factors for gluons would change, in such a way that Feynman diagrams for on-shell amplitudes 
end up being the same. 

To use the interaction with a background classical field as a matching condition, we must 
compute the current J M that couples to the background field A^ in Eq. (4. 10). Current conservation 
requires 

k-J(k) = Q, (4.11) 

where k is the external gluon's momentum. The usual convention for A^(k), from Eqs. (4.8) 
and (4.9), yields a current satisfying 

fc-J(fc) = 0, (4.12) 

where k \ = (2/a) sm(k^a/2). One sees, therefore, that a classical gluon line with Lorentz index jj, 
must be multiplied by 

One should think of n^k) as a wave-function factor for the external line. Its appearance has been 
noted previously by Weisz [17]. 

In the rest of this subsection we match the vertex function in lattice gauge theory with our new 
action to that in the continuum gauge theory. The incoming quark's momentum is p, the outgoing 
p', and the gluon's K = p' — p. The current is given by (no implied sum on fi) 

J, = n ll {K)N{p')u{i\p')k^p\p)u^p)N{p) ) (4.14) 

where A^(p',p) is the vertex function derived in Appendix A. The external quarks take normal- 
ization factors Af as well as spinor factors [7]. 
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1. Chromoelectric field: Li = 4 



For the interaction with the chromoelectric background field, we use the time component J4. 
To 0(p 2 /m 2 ) the current in continuum QCD is 

K 2 - 2zS • (K x P) 



J 4 = «(£', 0) 



8m 2 



w(£,o), 



where P = (jpf + f>)/2. After a short calculation with the new lattice action we find 

K 2 - 2iS • (K x P) z E K 2 a 2 ' 



J 4 = u(£',0) 



where 



1 



Am 2 E a 2 



8m| 



+ 



mo a 



u(£,0) 



C 2 ce 



+ 



2r, 



(4.15) 



(4.16) 



(4.17) 



[m a(2 + m a)] 2 m a(2 + m a) 1 + m a 
The correct (tree-level) matching is achieved if one adjusts 

z E = (4.18) 

and (ce, r E ) such that m E = m 2 : 

3 2m a(2 + m a) C 2 (C 2 — 1) ^,C 3 . r 2 ( 2 m a(2 + m a) 

(c E + r E — = t^— V + T^ + 77H \2 • (4 - 19) 

1 + m a m a(2 + mod) 1 + m a 4(1 + moa) 

At fixed m a the latter prescribes a line in the (c E , rE) plane. As before, this line becomes vertical 
at m a = 0, fixing c E = 1 and leaving r E undetermined. 

To obtain conditions on cee, tee, and Zee, we shall have to turn to Compton scattering in 
Sec. IV D. 



2. Chromomagnetic field: 11 = i 

For the interaction with the chromomagnetic background field, we use the spatial compo- 
nents Jj. To 0(p 3 /m 3 ) the current in continuum QCD is 



1 P 2 + \K 2 \ KiP-K 

m 



2m 3 ) 8m 3 



- s^Kj (± - P2 4 + J K2 ) + ,v/>:,/' ; ^} 0). (4.20) 

After another short calculation we find 

t -ffrtlpf 1 p2 + \ R2 \ KjP-K ZEa'K.P-K 

J, = -IU(£ ,0)<Pi[- - 2 + 



m 2 2m| / 8m2m E 7712(1 + mod) 

+ ±™ Bl a 3 [P 4 K 2 - iQ P • K] - ^wb^ s^K^K 2 

- IwB^e^KjPii?: ■ P + \w x a 3 Xi 

- Iw^P^Pf + \K 2 ) + ^w'^e^K^K 2 + K 2 ) 

+ + ti/Ja^tEj^KSif + iff?) + (3Pf + iffj)] 

- £ ^ (2^ ~ + £ ^ P ^}<^ 0) ' (4 - 21) 
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where m 2 , m\, W4, and m% have been introduced already, and 



1 + (c *- r * K (4.22) 



m B a m<2.a 1 + m^a 
1 _ 1 r s (r s -c B )( 2 | 8(z 6 -z 7 )+4(r 7 -4) 



m B ,a 3 mfa 3 (1 + m a) 2 1 + m a 

= 4(r g -c B )C 3 (l + m a) + 16(c 2 - c 3 )C + 8r 7 

[moa(2 + moa)] 2 mo£t(2 + moa) 1 + moa' 

w B2 = w B3 + — - — , (4.25) 

m a(2 + m a) 1 + moa 

w Bl = w Bo - gftl ~ 4) , (4.26) 
1 + m a 

Cj?C ~ 4(c 5 - r 5 ) 

w B = r- , (4.27) 

1 + m a 

r g C-24c 4 + 16(2c 5 + r 5 ) 

w 4 = — — r . (4.28) 

4(1 + m a) 

The term w x o?X is discussed below. 

Comparing Eqs. (4.20) and (4.21), one sees that the first four terms match the continuum if 
m 2 = m 4 = vtle = m. The other terms do not match unless one adjusts c B = r s [7] and z B = 
[as in Eq. (4.18)] and, furthermore, demands = w' 4: = w Bl = w B2 = w B , 3 = w' B = 0: 

r 7 m a(2 + m a) 

c 3 = c 2 + — — , (4.29) 

C 2(1 + m a) 

r' 7 m a(2 + m a) 

z 3 - -j nn . z~, (4.JU) 

C 2(1 + m a) 

c 4 = ^r s C + |c B C + 2r 5 , (4.31) 

c 5 = \c B ( + r 5 , (4.32) 

27 = 2 6 + !(r 7 -r 7 ), (4.33) 

4 = r 7 . (4.34) 

Taken with Eqs. (4.6) and (4.7), these tuning conditions put eight constraints on the nine (non- 
redundant) couplings for interactions made solely out of spatial derivatives (and, hence, chromo- 
magnetic fields). To eliminate z 6 from the right-hand side of Eq. (4.33), and to obtain conditions 
on r BB and z BB , we shall have to turn to Compton scattering in Sec. IV D. 

Equations (4.29)-(4.34) make concrete several abstract features of Sec. II. If one would like to 
take c 4 to be redundant in Eq. (4.7), then one cannot take r 5 to be redundant here, and similarly 
for zq and r 7 or r' 7 . Also, a mistuned C5 — r 5 leads to w' B 7^ and a spin-dependent contribution 
[1 + \w' B m2a[Kl + Kj)a 2 ]eijiiEiKj/2rri2. The mismatch here is suppressed by A 2 in the HQET 
counting — as expected from Table II — and by a 3 in the usual Symanzik counting. 

The only undesired term in Eq. (4.21) not yet discussed is iwja 3 ^, where 

X = (zS x K) P 2 - (iS x P) P ■ K - P [zS • (K x P)] + (K x P) iS ■ P, (4.35) 

= 4r s ( 3 (l + m a) | 16c 2 C 
[m a(2 + m a)] 2 m a(2 + m a) 



19 



One cannot tune wx = 0. Fortunately, however, X = 0. A simple geometric proof is as follows: 
if, by chance, P is parallel to K, then setting P oc K one sees that the last two terms on the 
right-hand side of Eq. (4.35) vanish and the first two cancel. In the general case that P is not 
parallel to K, then K, P, and K x P are three linearly independent vectors. But one easily sees 
that 

K ■ X = P ■ X = (K x P) ■ X = 0; (4.37) 

thus, X = 0. Such identities are very useful in simplifying expressions for the Compton scattering 
amplitude. 

C. Quark-quark scattering 

To match the four-quark action, S qqqq , one must work out the quark-quark scattering amplitude. 
With the current J u derived in the previous subsection, this is a relatively simple task. The main 
new ingredient is the improved gluon propagator. For k 2 a 2 <C 1, one finds [17] 

D^{k) = nii {k)D c p\k)n v {k) [l + xa 2 k 2 } + 0(a 4 ), (4.38) 

where x is the redundant coupling of the pure-gauge action, cf. Appendix C and Ref. [19]. This 
approximation suffices for evaluating t-channel gluon exchange. Once the bilinear action has been 
matched correctly, the lattice amplitude (using, say, Feynman gauge) is clearly merely 

A* (12 -> 12) = Aont(12 -> 12) + xa 2 t a .h ■ U\ (4.39) 

where 1 and 2 label the scattered quark flavors, and both t a have uncontracted color indices. We 
find, therefore, that the tree-level couplings of S qqqq are, at most, proportional to x. They can be 
eliminated, at the tree level, by setting x — 0, with the added benefit of simplifying the gauge 
action Spa F 2 . 

Note, however, that the approximation in Eq. (4.38) and, thus, Eq. (4.39), breaks down for s- 
channel annihilation of heavy quarks. As discussed in Sec. HE, these interactions are suppressed 
for other reasons, so the four-quark operators needed to correct them may be neglected. 

D. Compton scattering 

The matching of Sees. IV A-IVC leaves four non-redundant couplings of the new action un- 
determined: zq, cee, zee, and zbb- To find four more matching conditions, we turn to Compton 
scattering. We shall proceed with the gauge-action redundant coupling x = 0. 

The amplitude is 

<t(qg - qg) = ^ZWMWMfaMn^k), (4.40) 

where e„ and e M are continuum polarization vectors, and M. a ^ v denotes the sum of Feynman dia- 
grams shown in Fig. 2. The factors n u (k') and n^k) appear in Eq. (4.40) to account for lattice 
gluons. With them one can verify that 

^2 e ^k)n^{k)n u {k)e v {k) = -D^k), (4.41) 

pol. 
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as usual. We find it convenient to associate these factors with the diagrams and introduce Mf^ v = 
n u {k')Mf u n^{k). Then 

M* = t fe t a n,(A ; 0Ar(p0«(e , ,p / )A,(p / ,g)5(g)A M (g,p) M (e,p)Ar(p)^(fc) 
+ ^^ 6 n,(A ; OW)w(e^pOA M (p / ,g / )^(g / )A At (g / ,p)w(e,p)Ar(p)r^ At (A;) 
- i{^^^ b }n,(A ; OAr(pOM(e^pO«^(p,fc,-A ; , )M(e,p)Ar(p)r^ At (A;) (4.42) 

+ t c V^(k, -k', -K)D CTp (K)n,(A ; / )W)n(r,p / )Ap(p / ,p)u(e,p)Ar(p)^(fc) 

where q = p + k = p' + k' ', q' = p — k' = p' — k, K = k — k' = p' — p. The propagator 5(g) and 
vertex factors A M , X^ u and Y^ u are defined in Appendix A. The gluon propagator, to the accuracy 
needed, is given in Eq. (4.38), and to the same accuracy the triple-gluon vertex is (with x — 0) 

v^k-k'-K) = zr fe >^(£0M^> CT wr 1 { 

<W(* + k ')^ 1 - T2^K 2 a 2 ) + ±K a (k* - k'^a 2 ] 

- 5 va [{k' - K)^(l - ±5^k 2 a 2 ) + ±k,{k' 2 - Kl)a 2 } 

- 6 aiM [(K + £0,(1 - ^<L,fc'V) - ±k' u (K 2 - k 2 a )a 2 ] } . (4.43) 

Note that the factors n a (K), etc., arise naturally. Note also that K ■ J = k ■ e = k' ■ e' = k 2 = 
k' 2 = 0, so most of the lattice artifacts in the vertex drop out. The remaining one is necessary to 
cancel a similar lattice artifact from the other diagrams, cf. Eqs. (BIO) and (Bll). 

We may choose the polarization vectors such that e' A = e 4 = 0. Then we need only focus on 
■M-mn- We have verified that .M44 is improved by (a subset of) the improvement conditions needed 
for A(qg — > qg) calculated with these polarization vectors. 

The present the results, let us introduce some notation. Write the momenta as 

P = (p> +p)/2, {AAA) 
R = {k + k')/2, (4.45) 
K = p' -p = k-k\ (4.46) 

so q = P + R and q' = P — R. Note that P = —iP± = 2mi + • ■ ■ is larger than the other 
momenta, and K = —iK^ = (p' 2 — p 2 )/2m 2 is smaller. Next separate the diagrams according to 
a color decomposition, 

M% = \{t\ t b }M, v + \[t\ t b ]Af, u , (4.47) 



k,a,[i p' k,a : fi p' k,a,fj, k',b 7 v k : a,fi k'.b.u 




FIG. 2: Feynman diagrams for Compton scattering in lattice gauge theory. 
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where the second term would be absent in an Abelian gauge theory. Finally, write 

3 n 



M ^ = E E Ro-^M^- 1 - 2 ^, (4.48) 



n=0 s=0 

and similarly for N^, where the superscript (n, r) denotes the power in 1/m and R . 

Most of these terms are well-matched with Eqs. (4.18), (4.19), (4.29)-(4.34). New matching 
conditions come from Mmn , Nmn , Mmn , and Nmn^ ■ The (n, r) = (3, 2) amplitudes are 

• 1(3,2) _ S mn 2g 3 Z EE S mn 

Am EE 1 + m a 

(3j2) ErrmA (4 5Q) 

^' n EE 



where 



8[C + |cijCmoa(2 + m a)] 2 4( 2 



m EE a 3 [m a(2 + m a)} 3 [m a(2 + m a)] 2 

!6c EE ( 8(c EE ( + r EE ) 



(4.51) 



m a(2 + m a)(l + m a) 1 + m a 

To match to continuum QCD one requires 

zee = (4.52) 

and the adjustment of (c EE , r EE ) so that m EE = m 2 . As with, say, (c E , r E ), at fixed m a the latter 
prescribes a line in the (c EE , r EE ) plane, which becomes vertical at m a = 0, fixing c EE = — ~ 
and leaving r EE undetermined. 

The (n, r) = (3, 0) amplitudes are 

2a 3 

= M§£> | matchcd - —{z BB + z 6 + r 7 - r BB - z' 7 )M mn , (4.53) 

M mn = 5 mn (R 2 - \K 2 ) - (R m - \K m )(R n + \K n ), (4.54) 

2a 3 

A^ 0) = ^ | matchcd - — (^6 + r 7 - r BB - z' 7 )N mn , (4.55) 

N mn = ZmnriRri^ ' R ^K r iYj ■ K) — ^{iYi n S mrs + zS m £ nrs ) R r K s , (4.56) 

where "matched" denotes terms (spelled out in Appendix B) that already match, if the conditions 
derived so far are applied. Equations (4.53) and (4.55) yield the new conditions 

z BB + z 6 - z' 7 = r BB - r 7 , (4.57) 
z 6 - z' 7 = r BB - r 7 . (4.58) 

Solving these, and noting z' 7 = r' 7 [Eq. (4.34)], we find 

z BB = 0, (4.59) 
z 6 = r BB +r 7 - r 7 , (4.60) 

which completes the set of conditions needed to match the new lattice action. 
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E. Matching Summary 



Equations (4.6), (4.7), (4.31)-(4.34), (4.59), and (4.60) can now be combined to yield 

m a(2 + m a) 



6Cci 
16Cc 2 



~C + (cbC + 6r 5 , 

4C 4 (C 2 ~ 1) 
[m a(2 + m a)] 2 



1 + m a 

C 3 [2C + 4r s (l + m a) - 6r s ( 2 /(l + m a)] 



2/-4 



+ 



3rX 



c 3 

c 4 
c 5 

^3 
^6 

4 



(1 + 

C 2 + 



+ 



moa(2 + moa) 
m a) 2 ' 2(1 + m a) 
r 7 m a(2 + m a) (r s - 



m a(2 + m a) 
8(r B B + r' 7 - r 7 ) + 
c s )C 2 (l +m a) 



(4.61) 
(4.62) 



3^3 



1/-2 



1 + m a) 2 



m a 



C 2(1 + m a) 

\c B C + r 5 , 

r' 7 m a(2 + mod) 

J 2(1 + m a) ' 



4m a(2 + m a) 



r7, 



7): 



r 



7- 



0. 



(4.63) 

(4.64) 
(4.65) 

(4.66) 

(4.67) 
(4.68) 
(4.69) 
(4.70) 



To run a numerical simulation, we would like to have as few new couplings as possible. The 
matching calculations verified the presence of several redundant directions. We may, therefore, 
take 



r 5 = r 7 

to all orders in perturbation theory. Hence 



r B B = 



Cl 

c 2 = c 3 



■-C 



c B - 



c 3 (c 2 - 1; 



rriQa(2 + m^a) 
6(1 + m a) 

C 2 [C 



(4.71) 

(4.72) 
(4.73) 



2r s (l + m a) - 3r s ( 2 /{l + m a)} 



[2m a(2 + m a)}' 



8m a(2 + m a) 



+ 



3r 2 C 3 



C4 
C 5 



16(1 + mod) 2 
\cbC, 



m a(2 + m a)r 2 ( 
32(1 +m a) 2 



- 1 



m a 



and 



^3 — ^6 — ^7 



From the chromoelectric interactions we require 



zbb = 0. 
77^2 and 



+ 



r 2 m a(2 + m a) 



m 2 , whence 

te 2m a(2 + m a) 



moa(2 + moa) 1 + m§a 4(l+moa) 2 C 2 1 + moa 



(4.74) 

(4.75) 
(4.76) 



(4.77) 



(4.78) 
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c BE 2 + m a(2 + m a) = -= — — ^- + — — r 

|moa(2 + moa)J z moa(2 + moa) 

C(r s C-l-m a) , 2 x 2 

+ 2m a(2 + m a) + ^ C +2^-^(1 + ^) 

+ r ^ m ° a(2 + m ° a) - ^m a(2 + moa), (4.79) 
1 + m a C 



and we also find 

z E = zee = 0. (4.80) 

Without loss one may set the redundant r E = t E e = to simplify the action and Eqs. (4.78) 
and (4.79). 

In summary, of the nineteen new couplings in Eqs. (3.5)-(3.7), we find only six that are non- 
zero at tree-level matching. Moreover, once the bilinear action has been matched, and the redun- 
dant gauge coupling x — 0, the only non-zero four-quark interaction would correspond to (highly 
suppressed) QQ annihilation. In the next section we shall examine the size of the remaining un- 
certainties, to justify that this level of matching suffices. 



V. ERRORS FROM TRUNCATION 

In this section we give a semi-quantitative analysis of heavy-quark discretization effects with 
the new action. Our aim is to study the accuracy needed in matching lattice gauge theory to 
continuum QCD. Several elements are needed. First, we need estimates of the mismatch at short 
distances. This is straightforward, because the calculations of Sec. IV can be applied to work out 
how large the mismatch is for the unimproved action. Second, we need estimates of the long- 
distance effects, which is possible parametrically, by counting powers of A and v. Finally, the size 
of discretization effects depends on the lattice spacing (obviously) so we must note the range that 
is tractable today and in the near future. 

The error analysis is convenient using the non-relativistic description. Heavy-quark effects of 
operators that are related as in Eqs. (2.14) and (2.15) are lumped into one short-distance coeffi- 
cient C\ at per HQET operator in Table III. In Sec. IV the short-distance coefficients are l/2m 2 , 
l/2m B , l/4m%, l/8m|, 104, w Bi , etc. In the corresponding continuum short-distance coeffi- 
cients Q ont , these masses are replaced with a single mass rriQ. To eliminate discretization effects 
from the kinetic energy, one should identify txiq with m 2 . 

Comparison of Eqs. (2.5) and (2.6) then says that heavy-quark discretization effects take the 
form 

error, = ($* - C?*) (O*) . (5.1) 
For example, the error from (p 2 ) 2 /&m\ is 

See Refs. [11, 12] for further details, and Ref. [31] for the application of this technique to compare 
several heavy-quark formalisms. We estimate the matrix elements (Oi) using the power counting 
of HQET and NRQCD for heavy-light hadrons and quarkonium, respectively. The power of A 
or v is listed in Table III. The coefficient mismatches are obtained from Sec. IV, where explicit 
expressions show how the coefficients depend on the new couplings. In particular, when the new 
couplings vanish, we derive the mismatch for the Wilson and clover actions. 
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Explicit calculations of the mismatch at higher orders of perturbation theory are not yet avail- 
able. (They would be tantamount to higher-loop matching.) Nevertheless, the asymptotic behav- 
ior remains constrained, when mqa <C 1 because of the presence of the £-type operators, when 
m,Qa 1 by heavy-quark symmetry, and when mqa ~ 1 because the Wilson time derivative 
ensures only one pole in the propagator [11]. It turns out that the most pessimistic asymptotic be- 
havior for l/2m B , l/Am 2 E , etc., is the same at higher orders as in the tree level formulas in Sec. IV. 
It seems reasonable, therefore, to multiply the tree-level mismatch with a l s to estimate the /-loop 
mismatch. We use one-loop running for a s (a) starting with a s (l/ll fm) = 1/3. This yields the 
high end of the Brodsky-Lepage-Mackenzie coupling [32] calculated for similar quantities [33]. 

The resulting estimates for the mismatch of rotationally symmetric operators are shown in 
Fig. 3, as a function of the lattice spacing a = m2a/m,Q, Q G {c, b}. We show the relative error 
in mass splittings, which are of order A in heavy-light hadrons and of order mqv 2 in quarkonium. 
The left set of plots uses HQET power counting, for heavy-light hadrons, while the right set of 
plots uses NRQCD power counting, for quarkonia. The light gray or red (dark gray or blue) curves 

HQET for heavy-light NRQCD for quarkonia 



relative error relative error 




a (fm) a (fm) 

FIG. 3: Relative truncation errors for the new action. The light gray or red curves stand for c quarks; 
dark gray or blue for b. Dotted curves show the error when the contribution is unimproved. Dashed and 
solid curves show the error for tree-level and one-loop matching, respectively, of the needed operators. 
A = 700 MeV, m c = 1400 MeV, m b = 4200 MeV; v\ c = 0.3, v\ b = 0.1. Vertical lines show lattice 
spacings available with the MILC ensembles [34]. 
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show the estimate for hadrons containing c (b) quarks. The dotted curves show the error when the 
corresponding correction term is omitted completely, i.e., the errors in the Wilson action. The 
dashed (solid) curves show the estimate of the error for tree-level (one-loop) matching. The verti- 
cal lines highlight a = 0.125 fm, 0.09 fm, 0.06 fm and 0.045 fm, corresponding to the ensembles 
of gauge fields with n j = 2 + 1 flavors from the MILC collaboration [34] . 

To drive the each contribution to heavy-quark discretization effects below 1%, we find that one- 
loop matching is necessary for cb, the coupling of the chromomagnetic clover term. Tree-level 
matching is sufficient for the chromoelectric clover coupling ce, though one-loop matching would 
be desirable for charmonium and charmed hadrons. The lowest plots, labeled "from l/Smf" are 
for the relativistic correction terms, with couplings C2 and z 6 . They also apply to l/8m B , and the 
related chromomagnetic couplings C3 and z 7 . The one-loop mismatches of four-quark interactions 
are suppressed not only by a loop factor, but also by A 2 or v 2 , so they should fall below 1% too. 

Similar results for operators that break rotational symmetry are shown in Fig. 4. To drive these 
contributions to heavy-quark discretization effects below 1%, we again find it sufficient to tune the 
couplings of the new action at the tree level. 

There are some other noteworthy features of Figs. 3 and 4. For rriqa <C 1, the discretization 
effects vanish as a power of a, as one would deduce from the Symanzik effective field theory. 
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NRQCD for quarkonia 
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FIG. 4: Relative truncation errors for the new action, from discretization effects that break rotational sym- 
metry. The curves have the same meaning as in Fig. 3. 



26 



Because we identify m 2 with the mass in the Cf ont , the powers of a are balanced by A or rriQV, 
not rriQ. Had we identified mi with the physical mass, errors of order {rriQa) n would have ap- 
peared. For rriQa ~ 1, the tree-level curves flatten out. The error cannot grow without bound, 
because of the heavy-quark symmetries of the Wilson action and our improvements to it. Indeed, 
the curves for the b quark are usually lower than those for the c quark, which bodes well for calcu- 
lations relevant to the Cabibbo-Kobayashi-Maskawa matrix. The underlying reason for the pattern 
is that the static approximation works better for ^-flavored hadrons than for charmed hadrons. The 
1/m^ contributions start out smaller, so their mismatches are also smaller. Similarly, the leading 
NRQCD works better for bottomonium than charmonium. The mismatches from 1 / 8m| and W4/6 
deviate from the pattern, however, because NRQCD's relative suppression v^Jv%, is not as strong 
as HQET's (m c /mb) 3 . Mismatches from wbJ^ and (w^ + w'^/A are of order v 4 and again follow 
the pattern. 

In tree-level improvement, one should avoid choices where it is known that one-loop correc- 
tions from tadpole diagrams will be large [35]. Therefore, we envision following some sort of 
tadpole improvement. In the action, write each link matrix as uqIU^/uq] and absorb all but one 
pre-factor of u into a tadpole-improved coupling Cj and f (In several cases, it will be neces- 
sary to expand expressions such as D iut A ilat , A^ at , and Eq. (3.11), to eliminate any instance of 
UpU^ = 1 before inserting u .) Then apply the conditions of Sec. IV to q and fj instead of q 
and 7"i, and take the uq factors in the denominator from the Monte Carlo simulation. 

VI. CONCLUSIONS 

In this paper we have presented the formalism and explicit calculations needed to define a new 
lattice action for heavy quarks. Our aim was to obtain an action whose discretization errors would 
be < 1% at currently available lattice spacings. Combining our matching calculations, power 
counting, and the heavy-quark theory of discretization effects, we have argued that the proposed 
action should meet its target. Setting to zero the redundant couplings and those that vanish when 
matched at the tree level, our action can be written S = S + S B + S E + S'ncw* where 

S^w = c 1 a 6 ^2ij(x)^2-f l D ilat A il£Lt ip(x) +c 2 a 6 ^^(x){7- D lat , A\H}ip(x) 

x i x 

+ C 3 a 6 ^^(x){7 • Aat,«S ■ B lat }?p(x) + C££;a 6 ^^ (x){7 4 -D41at,« • Ei at }ip(x) 

X x 
x i x i jj^i 

The new action has six additional nonzero couplings, which depend on the couplings in Sq + Sb + 
Se according to Eqs. (4.73)-(4.76) and (4.79). To achieve 1% accuracy, Sb must be, and Se could 
well be, matched at the one-loop level [36]. 

Another lattice action achieves similar accuracy for charmed quarks, namely the highly- 
improved staggered quark (HISQ) action [37]. Our approach is computationally more demanding 
than HISQ. Its advantage, however, is the intriguing result that our discretization errors for bottom 
quarks are smaller than for charmed quarks. That means that experience with charmed hadrons 
and charmonium can inform analogous calculation of properties of b- flavored hadrons. 

Finally, we note that there is tension between the most accurate calculation of the D s meson 
decay constant, fu s [38], which uses HISQ, and experimental measurements [39]. Our action is a 
candidate for the charmed quark in a cross-check of the HISQ fo s , because its discretization errors 
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can be expected to be small enough to strengthen or dissipate the disagreement, while possessing 
different systematic errors. 
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APPENDIX A: FEYNMAN RULES 

In this Appendix we present Feynman rules for the new action needed to carry out the matching 
calculations of Sec. IV. These are the quark and gluon propagators and three- and four-point 
vertices. The corresponding Feynman diagrams are shown in Fig. 5. 

The quark propagator [Fig. 5(a)] is modified only through c 2 , c\, z 6 , and C4. It reads 

aS' 1 (p) = ry4 sin(j»4a) + 27 • K(p) + fx(p) — cosfj^a) (Al) 

where 

Ki{p) = sm{ Pi a) [C - 2c 2 pV - c x p 2 a 2 \ (A2) 

Hip) = 1 + m a + pV [|r s C + z 6 p 2 a 2 ] + c 4 ^fea) 4 (A3) 

i 

The tree-level mass shell is p 4 = iE, where the energy satisfies 

cosh Ea = 1+ ^ 2 + K2 . (A4) 

Mp) 

Incoming external fermion lines receive factors u(£, p)M{p) or v(£, p)J\f(p), where 

"<») = (^yiin) 1 ' 2 - < A5 > 

/>. \ L + sinh E — i~f ■ K .. . 
L + sinh E + ry • K , 

= TfTW*'"' (A?) 

L = fi(p) — cosh£J; 74«(^, 0) = u(£, 0), 74u(^, 0) = —v(£, 0). Outgoing external fermion lines 
receive factors M (p)u(£,p) or Af(p)v(£,p), where tt(£,p) = w^(^,p)74 5 v(£,p) = v^^p)^^. 

The gluon propagator [Fig. 5(b)] is not easy to express in closed form. We refer the reader to 
two papers of Weisz for details [17] and a correction [18] for the propagator on isotropic lattices. 
The improved vertex is in Ref. [18]. 
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Now let us turn to vertices with one [Fig. 5(c)-(d)] or two [Fig. 5(e)-(g)] gluons attached to a 
quark line. The new terms in the bilinear part of the action are all built from difference and clover 
operators that already appear in So + Sb + Se- Consequently, the new terms in the Feynman rules 
for these vertices can be obtained using the chain rule. 

The difference operators are given in Eqs. (3.8)— (3. 10). To simplify notation, let us drop the 





FIG. 5: Feynman rules for the action S given by Eqs. (3.1)— (3.7). 
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subscript "lat" in this Appendix. One-gluon vertices need 



Ap*(P,A) 



3D, 



dAlik) 



g t a S PPi cos[(P + \k) p a b 



dA, 



dA'Jk) 



g t a S m (2i/a) sm[{P + \k) 



dF 

p a (k)= pa 
dA'Jk) 



g t a cos \k p a [5 pa iS p (k) - 5 pp iS a (k)] . 



It is convenient to write out the chromomagnetic and chromoelectric cases of Eq. (A 10): 



Bi a m (k) 



-g t a cos(±k m a)e mri iS r (k), 



dA« m (k) 
dEi 

Ei a m (k) = — * M = O £ a cos(§A; m a)<S mi z&i(A;). 



BE; 



—got a cos(~k4,a)iSi(k) , 



dAi(k) 

since Bj = \eijkFj k and i?, = F 4i appear in Eq. (3.1). Two-gluon vertices need 

d 2 D n 



(A8) 
(A9) 
(A10) 

(All) 
(A12) 
(A13) 



A,£(P,M) 



«9 2 A 



rp ab 



(k,l) 



d 2 F 



pa 



dA*Jk)dA*(l) 



Ei,t(k,l) 
Et£(k,l) 



d 2 B; 



dA« m (k)dA» n (l) 

d 2 E t 
dA" m (k)dA* n (l) 

d 2 E t 
' dAl{k)dA b n {l) 
d 2 E t 



(A14) 
(A15) 

(A16) 



9l[t a A {{SwSu* - 5 p J vp )C pv {k,l) (A17) 
\5 pv a 2 K p [6 PP (S a (k) - S a (l)) - <V (S p (k) - S p (l))]} , 

9o\P ' ^ ] | ^ vaniC m rSk i '^mn^mri®' K m \S r (Ji) S r (l 



where K = k + /. For the clover operator it is convenient to introduce 

Cfiu(k, I) = 2 cos \{k + l) p a cos \l p a cos \{k + t) v a cos \k v a — cos \k p a cos ~J> v a. 
Then one has (K = k + /) 



% 2 [*V 6 ]l<W mi a 2 £ m [^(fc) - 5 4 



dA%{k)dA\{l) 
The Feynman rules for one gluon are then 

Fig.5(c,d) = -so^A>>) 
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(A18) 
(A19) 

(A20) 

(A21) 

(A22) 



with 



A 4 (p',p) — 7 4 cos[|(p' +p) 4 a] - isin[~(p' + p) 4 a] + |c£Caa ■ S(k) cos(|/c 4 a) 
+ zr B a 2 74 S ■ {S(k) x [S(p') + S(p)}} cos(|fc 4 a) 

- (r s - z E )a 2 'j A S(k) ■ [S(p') - S(p)} cos(|/c 4 a) 

+ ceeo 2 ! ■ S(k) [5*4 (p') - 5* 4 (p)] cos(|A; 4 a), (A23) 

Kn(p',p) = C^ m cos[\(p' + p) m a}- ir s (s\n[\(p' + p) m a\ 

- ^CB(ae mr iT>iS r (k) cos(|A; m a) — IceCaamS^/c) cos(|/c m a) 

- irE^emriEa^S^k) [S r (p') + S r (p)} cos(\k m a) 
+ (r E - z E )a 2 ^S^{k) [^(p') - S m (p)] cos{\k m a) 

- c 2 a 2 { 7m cos[|(p' + p) m a] (p' 2 + p 2 ) + 7 ■ [S(p') + S(p)] (p'Tp)™} 

- |cia 2 7 m |cos[|(p' + p) m a] (p^ 2 + p^J + [S m (p') + S m (p)] (p' + p) m } 

- c 3 a 2 £ m rt7475»SV(fc) [<Si(p') + cos(ifc m a) 

+ (c 3 - z 3 )a 2 7 ■ S(k) [S m (p') - 5 m (p)] cos(|fc m a) 

- (c 3 - 2 3 )a 2 7 m 5(/c) • [S(p') - 5(p)] cos(|fc m a) 

- CEEa 2 lmS±{k) [S 4 (p') - S^p)} cos(±k m a) 

- iz 6 a 3 (p r +p) m (p /2 +p 2 ) 

- ^c 4 a 3 (p'+p) m (p'J +p 2 m "j 

- (z 7 + c 5 )a 3 e mri T,iS r (k) (p /2 + p 2 j cos(|/c m a) 

+ c 5 a 3 £ m HSi5' r .(A;) (p'i + p 2 ) cos(|/c m a) 

+ r 5 a 3 £ mri £i,SV(fc) [Si(p')Si(p)] cos(\k m a) 

+ (r 7 - 4 - r 5 )a 3 e mri SiS' r (/i;) [S(p') • 5(p)] cos(±fc m a) 

- r 7 a 3 e mri [Si(p')H • 5(p) + S;(p)£ • 5(p)] SV(A;) cos(|A; m a) 

+ ,( r7 _ r ^) a 3 [5 m (p')5(p) • S(k) - S m (p)S(p') ■ S(k)} cos(±k m a). (A24) 

In the r 5 and z' 7 terms, Eq. (3.1 1) has been assumed. If instead one prefers Eq. (3.12) then replace 

[SjWM - [cos(|%a)p^] . 

Both choices have the same effect on Eq. (4.21). 
The two-gluon rules are 

Fig. 5(e, f , g) = -\gl{t\ t h } t3 aX, u {p, k, I) - \gl\t\ t h ] l3 aY, v {p, k, I), (A25) 

with 

X mn (p,k,l) = i(5 mn j m sin(y m a) - r s (5 mn cos(ls m a) 

- 2r E ae mn aJli [cos{\s n a) cos{\k n a) cos{~k m a)S±{k) 
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— cos(|s m a) cos(hl m a) cos(|/„a)5 , 4(/)] 

+ i(r E - z E )a 2 -f 4 SmnSm(^s m a) [Sm^S^k) + S , m (/)5 , 4 (/)] 
+ 4ic 2 7 m [cos(±s m a) cos(|Z m a) sin(±s n a) cos(|/c n a) 

+ sin(|s. m a) sin(|/ m a) cos(^s n a) sin(i/c n a)] 
+ 4ic 2 7n [sin(~s m a) cos(|Z m a) cos(~s n a) cos(|fc n a) 

+ cos(|s m a) sin(|/ m a) sin(is n a) sin(|A;„a)] 
+ 2ic 2 a5 mn cos(\s m a) 7 ■ [S(p') + S(p)] 

- ic 2 a 2 (5 mn 7 m sm(|s m a) (p /2 +P 2 ) 

+ icia5 mn -i m s m [4cos(|s m a)cos(|fc m a)cos(|/ m a) - l] 
+ 2ic 3 e mnr ^ 5 [sin(i r a)cos(|s TO a)cos(|i m a)cos(|i n a) 

— sin(A; r a) cos(^s n a) cos{\k n a) cos(|fc m a)] 

+ 2z(c 3 - z 3 )a {[<5 m „7 • 5(7) - 7„S , m (Z)] sin(|s m a) sin(|/ m a) cos(|/ n a) 
+ [<5 m „7 ■ S(k) - 7 m S n {k)\ sin(|s n a) sm(\k n a) cos(|/c m a)} 

- 8z 6 [sin(|s m a) cos(|i m a) sin(|s„a) cos(|A; n a) 

— cos(|s m a) sin(|/ m a) cos(|s n a) sin(|/c n a)] 

- 2z§a 2 5 mn cos(\s rn a) (p' 2 + p 2 ^j 

- 2c A a 2 5 mn |cos(|s m a) (p' 2 m + p 2 7 J + cos[|(A; - l) m a]s 2 m - k m l m \ 

+ 2i(z 7 + c 5 )a 2 Si [s n £ mr i5 r (A;) cos(i/c m a) cos(ifc ra a) 

+ s m e nri 5 r (0 cos(|/ n a) cos(|/ m a)] 
+ 2zc 5 a 2 £ mm . [s n S n 5 r (A;) cos(|A; m a) cos(|A; n a) 

— SmE m S r (l) cos(~Z n a) cos(|/ m a)] 

+ ir 5 a 2 e mnr {T, n S n (s)S r (k) - E m S m (s)S r (l)} cos(|fc m a) cos(|/ n a) 
+ 2r 7 a 2 S n e mri S , r (/c) {[<%(p') + Si(p)] cos(±s n a) cos(±fc n a) 

+ [Si(p') - Si(p)] sin(is n a) sin(|/c n a)} cos(|A; m a) 
+ ir 7 a 2 H m e nri S r (l) {[^(p') + S^p)] cos(|s m a) cos(|Z m a) 

+ [<Si(p') - Si(p)] sin(|s m a) sin(±/ m a) } cos(|/ n a) 

- ir 7 a 2 s mnr S r (k)Y, • {[S(p') + S(p)} cos(|s n a) cos(±£; n a) 

+ [S(p) - S(p)] sin(|s n a) siri^A^a)} cos(\k m a) 
+ zr 7 a 2 £ mm ,SV(/)£ • + S(p)] cos(|s m a) cos(|/ m a) 

+ [S(p') - S(p)]sin(|s m a) sin(|/ m a)} cos(i/ n a) 
+ z(4 + r 5 - r 7 )a 2 e mri S n (s)T,iSr(k) cos(|A; m a) cos(|/„a) 
+ i(z 7 + r 5 - r 7 )a 2 e nri S' m (s)S i S , r ,(/) cos(|/c m a) cos(±/ n a) 

- (r 7 - r 7 )a 2 S n (k) [S m (p') - S m (p)} cos(|s n a) cos(ifc n a) cos(|/c m a) 

- (r 7 - r 7 )a 2 S m {l) \S n (p') - S n (p)} cos(|s m a) cos(|/ m a) cos(~Z n a) 

- (r 7 - r 7 )a 2 S n (k) \S m (p') + 5 m (p)] sin(±s n a) sin(±fc n a) cos(\k m a) 
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- (r' 7 - r 7 )a 2 S m (l) [S n (p') + S n (p)} sin(±s m a) sin(±Z m a) cos(|Z n ,a) 



r 



- r 7 )a 2 5 mn S(k) ■ [S(p') - S(p)] cos{\s m a) cos 2 (\k m a) 



+ (r' 7 - r 7 )a 2 5 mn S(l) ■ [S(p') - S(p)} cos(|s m a) cos 2 (|/ m a) 

+ \{r' 7 - r 7 )a 4 6 mn s m {S m (k)S(k) + S m (l)S(l)} ■ [S(p') + S(p)] 

+ 2(r BB - z BB )a 2 [5 mn S(k) ■ S(l) - S m (l)S n (k)] cos(\k m a) cos(±Z n a) 

- 2(r EE + ZEE)a 2 5 mn S±{k)S±(l) cos(\k m a) cos(|Z n a), (A26) 

where now p' — p + k + Z, and s = p' + p = 2p + k + Z; 

X 44 (p,kJ) = i^ 4 sm[\(p' + p) 4 a] - cos[|(p' + p) 4 a] 

+ icEEO 2 [7 • S(k)S 4 (k) + 7 • S(l)S 4 (l)} sin[|(p' + p) 4 a] 

— 2{jee + ZEE)o 2 S(k) ■ S(l) cos(|/c 4 a) cos(|Z 4 a), (A27) 
X 4m (p,k,l) = -2r E ae mr aJliS r {k) cos(±s m a) cos(±/c 4 a) cos(|/c m a) 

— 2(r£ — ZE)a 2 7 4 /c^ sin(|s m a) cos(|/c 4 a) cos(|/c m a) 

— icEEa 2 lmit sin[~(p' + p) 4 a] cos(|Z 4 a) cos(±Z m a) 

+ 2(r E E + ZEE)a 2 S m (k)S 4 (l) cos(|/c 4 a) cos(|Z m a), (A28) 

^mn(P, fc, Z) = -iCsC^iCmni^k, Z) - |c E Ca 2 5 mn Q: m i^ m [S 4 (&) - S^Z)] 

- |r jB a 3 £ m „ i 7 4 E i s n k n S 4 (k) cos(~A; m a) + s m l m S 4 (l) cos(~Z n a) 

- |r E a 3 5 mn £: mr . i 7 4 S i K m [S' r .(p / ) + S , r (p)][5 , 4 (A;) - S^Z)] 

+ 2z(r£; - z E )a8 mn r y 4 [S 4 (k) cos 2 (|/c m a) - 5* 4 (/) cos 2 (|Z m a)] cos(|s m a) 

- |(r B - z E )a 3 j 4 5 mn K m [S m (p') - S , m (p)][5 , 4 (/c) - S^Z)] 

- 42c 2 7 m [cos(|s m a) cos(|Z m a) cos(|s n a) sin(|/c n a) 

+ sin(|s m a) sin(|/ m a) sin(|s n a) cos(|fc n a)] 
+ 4zc 2 7n [cos(|s n a) cos(|/c n a) cos(|s m a) sin(|Z m a) 
+ sin(|s n a) sind^a) sin(|s m a) cos(|/ m a)] 

- 2c 1 5 mn ij m cos[(p' + p) m a] sin[±(A; - Z) m a] 

- 2ic z a^ 4 ^C rnni {k,l)[S i {p') + ^(p)] 

4 c 3 a 7475^mnr krknS n COs( 2 k m a) + l r lmSm COS^Zraft) 

- 2z(c 3 - z 3 )a K7 mre (A;, Z) {7 m [S' n (p / ) - 5 n (p)] - 7„[5 m (p') - S m (p)]} 

+ ±5 mn a 2 iU>W) - 5 m (p)] 7 • [S(k) - S{1)} 

- \5 mnlm a 2 k m [S{p') - S{p)] ■ [S(k) - S(l)} 
+ [8 mv n ■ S(l) - j n S m (l)] cos(|s m a) cos(|Z m a) cos(±Z re a) 

- [5 m „7 ■ S(k) - jmSn(k)} cos(is n a) cos(±A; n a) cos(±/c m a)) 
+ |c£ S a 3 7 m (5 mn K m [S 4 (p') - S 4 (p)} [S 4 (k) - S 4 (l)} 

- 2z 6 a 2 cos(§s m a)Z m s n cos(|A; ri a) - cos(\s n a)k n s m cos(\l m a) 
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+ 4c 4 5 mn sin[(p' + p) m a] sm[\(k - /) m a] 

- 2z(z 7 + c 5 )a 2 £; e mri S r {k)k n cos{\s n a) cos(|/c m a) 

- £nriS r {l)lm COs(|s m a) COs(|i„a) + C mni (k, I) (j/ + p 



T, n S r (k)k n cos(|s n a) cos(|/c m a) 



+ E m 5V(0^mCOs(|s m a) cos(i/ n a) 
+ 2ic 5 a 2 Si (p'J + p 2 ) C mni {k, I) 



-\- ir^a £ mnr Yl n 



K n cos(|/ n ,a) — l n sin 2 (|s n a) S r (k) cos(\k m a) 
K m cos(±/c m a) - k m sin 2 (|s m a) S r (l) cos(±Z n a) 



+ 2ir 5 a 2 T li S i (p')S i (p)C rnni (k,l) 

+ ir 7 a 2 T, n e mri S r (k) {[Si(p) - S^p)} cos(|s„a) cos(|/c n a) 
+ + Si(p)] sin(|s n a) sin(±/c n a)} cos(|/c m a) 

- ir 7 a 2 Y, m e nri S r (l) {[<%(?/) - Si(p)] cos(is m a) cos(i/ m a) 

+ [<%(£>') + Si(p)] sin(|s m a) sin(|/ m a) } cos(|/ n a) 

- ir 7 a 2 e mnr S r {k) {£ • [S(p') - cos(~s m a) cos(|/c m a) 

+ X • [S(p) + S(p)] sin(±s m a) sin(|/c m a)} cos(|fc m a) 

- ir 7 a 2 e mnr S r (l) {S ■ - cos(|s n a) cos(|/ n c 

+ S ■ [£(p') + S(p)] sin(is n a) sin(±Z n a)} cos(|/ r 

- 2zr 7 a 2 £ • [S(?/)^(p) + S(p) C mm (k, I) 

+ i(z' 7 + r 5 - r 7 )a 2 e mri 5' r (A;)S i jif n cos(|/ n a) -/„ sin 2 (|s n a)j cos (±/c m a) 

- z'(4 + r 5 - r 7 )a 2 £ nri S r (l)Ei ^k m cos(\k m a) - k m sin 2 (|s m a)| cos(|Z n a) 

- 2i(z 7 + r 5 - r 7 )a 2 Si5(p') • S(p)C mni (k, I) 

- (r 7 - r 7 )a 2 [S m (p') + S m (p)] S n (k) cos(|s n a) cos(|fc re a) cos{\k m a) 
+ (r 7 - r 7 )a 2 \S n (p') + S n (p)} S m (T) cos(|s m a) cos(|/ m a) cos{\l n a) 

- K - r 7 )a 2 [S m (p') - S m (p)} S n (k) sin(±s n a) sin(|A; re a) cos(|A; m a) 
+ (r' 7 - r 7 )a 2 [S n (p) - S n (p)} S m (l) sin(±s m a) sin(|/ m a) cos(|/ n a) 

ri)a 2 5 rnn S(k) ■ \S(p') + S(p)] cos(±s n a) cos 2 (\k n a) 
r 7 )a 2 5 mn S(l) ■ [S(p') + S(p)} cos(±s m a) cos 2 (|/ m a) 



+ 



{S m (k)S(k) - S m (l)S(l)} ■ [S(p') - S(p)} 



+ 2(r' 7 - r 7 )a 2 [S m (p')S n (p) - S m (p)S n (p')} C mn (k, I) 



r 7 )a 4 5 mn K m [S m (p')S(p) - S m (p)S(j/)] ■ [S(k) - S(l)] 
+ irBB^^mnr [S r (k)H ■ S(l) + S r (l)H ■ S(k)] cos(|/c m a) cos(|Z n a) 
+ ir BB a 2 (S m e nri + S n e mri ) S r (k)Si(t) cos(|/c m a) cos(±Z n a) 
- 2ir E Ea 2 E mni Y H S±{k)S±{l) cos(\k m a) cos(±/ n a), 



(A29) 
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where C mn i(k, I) — e mn iC mn (k, I) ^5 mn e mr ia 2 K m [S r (k) S r (l)]', 

Y AA {p,k,l) = ~c E (,aa ■ [S(k) - S(l)] sin[|(A; + Z) 4 a] 

- r E a 2 lA H ■ {[S(p') + 5(p)] x [S(k) - S(l)}} sin[|(A; + /) 4 a] 
+ z(r s - ^)a 2 74 W) - S(p)] ■ [S(k) - S(l)} sin[±(A; + Z) 4 a] 

+ 2icEE a [7 ■ cos 2 (|/c 4 a) — j ■ S(l) cos 2 (|Z 4 a)] cos[|(p' + p) 4 a] 

- 2ic EE aj ■ [S(k) - S(Z)] sin 2 [±0 + Z) 4 a] cos[|(p' + p) 4 a] 

- 2zr Bi; a 2 £ ■ [S(Jfe) x 5(Z)] cos(±£; 4 a) cos(±Z 4 a), (A30) 
Y Am (p,k,l) = -c E (a m C4, m (k,l) 

- 2r E ae mr aJ^i[Sr(p) + S r (p)]C 4m (k,l) 

- r E a 2 £mHlJ^i sm{\s m a)k m S r {k) cos(|/c 4 a) 

- 2i(r E - z E )a^ A [S m {p') - S m (p)}C Am (k, I) 

- 2i{r E — z E )a^ A S m {k) cos(|s m a) cos{\k A a) cos(|Zc m a) 
+ 2ic EE a-f m [S 4 (p') - S 4 (p)]C 4m (k, I) 

+ 2ic EE a^ m S 4 {l) cos[|(p' + p) 4 a] cos(|Z 4 a) cos(|Z m a) 

- 2ir EE c?e mr iTiiS r {k)Si{l) cos(|/c 4 a) cos(|Z m a). (A31) 



APPENDIX B: DETAILS OF COMPTON AMPLITUDES 



The parts of the Compton scattering amplitude not exhibited in Sec. IV D are shown here. First 
the color-symmetric contributions: 

M$£> = 6 —, (Bl) 

• >(2 —1) -Pm(R ~\~ ~^ Rn(R 2~^) m 

mn ' = m\ 

[Cg - \K) m e nri (R - \K) r -(R + \K) n e mn (R + ^) r ]zE t 

2mim E 

2( y iYj rn £ TlTS -\- iYj n £ mrs ^R r K s A£ mnr R r i^l • R -\- £ mnr K r i^l • K , D », 

+ 8ml ' (B2) 

M ™ = (B3) 
= [AP m P n + (R- \K) m {R + \K) n ] m \~Jf* 

+ [P m (R+\K) n + P n {R-\K) m ]?—^- 

m 2 

AR 2 — K 2 

- [e mri (R + \K) r (R + \K) n - e nri (i? - |JOr(i? - \K) m ] 1 



'2m\mE 
AR 2 — K 2 

\_i.R 2^) m ("^ 2^") n ~2^\^-"n^ iTir s i^m^nrs)RrK- s\ ~^Q^~^2 
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(AR 2 — K 2 ) 2 P ■ R 

bAm,2m B 

~t~ \^rnnrK r t^j ' R, 8 mnr R r 1^j • _K" 



{^^wnnrRr^- 1 ' E-rnnrK r 7^j ' -K- 



M {m \ 

mn | match 



A P 2 , op p 



2m\ 



1 



1 



m\ m B m E 
2z E a 2 



' Am2m B 

AR 2 - K 2 

32m 2 m 2 3 
PR 
8m2'm 2 B ' 
1 \ 4i? 2 + K 2 

2~ j 



(B4) 



16 



+ a 



4m 2 m|; 4m B m 2 E e mia m2 
r 2 - 4)C 2 



i^RmRn ^K m K n ) 



> 2 - 4)c 2 



+ 



2m 3 3 , 2m 2 m|, 



+ a 3 ^ 2 J 5 mn (4# 2 -K 2 ) 

+ a 3 ^^ (4i? mJ R„ - K m K n ) 
(2P 2 + ij 

&m.rt.i^- J iP ' R^ 



+ a 3 lw Bl (SmnK 2 - K m K n ) - aWmn (2P m + |i? 2 , + J2 Km) 
1 1 



+ a 2(1114 + u£ 



5 =- + ja 3 ws 3 I e mnr P r iYl • R 

Am2m E Am B m E J 

_|_ _|_ 3 1 / _|_ / \ 3 3 

2m 3 3 , 4m2m E 4m B m 2 E 2 443 



^m.n.rRrT^ ' P 



— a 



[rl ~ 4)C 2 , 



gg2mia 



(RmK n R n K n 



8m2fflg 
1 



(K n £ mr i -\- K m E nr i) P T iYji 



2 O g^i?3 ) (^S n £ mrs -|- %Yj m £ nrs }P r K s 
oTYl B 7Tl E 



4m B , %m27v? E 



+ W(w 4 + w' 4 ) 



— a 4(1^4 + w ^£ mnr K r {^P m iYi m — P n iX n ). 
The color- antisymmetric contributions from Fig. 2(a)-(c): 



JV(1,0) 
* mn 

^(2,-1) 



4P m P n + (P-|if) m (P+|if) n 
2m 2 , 

[Pm^nri(P 2 ~^0 r Pn^mri(P "I - 2 -^Or]^^i 

m 2 m B 

4ml 



(B5) 



(B6) 
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6 mn (R 2 - \K 2 ) -(R- \K) m {R + \K\ 



2m% 



5 mn Aa 2 z E 8 



2mi 1 + m a 



(B7) 
(B8) 



PR 



^n~ 2) = ~ [4PmPn +(R~ \K) m {R + \K) n ] ^-3 



[P m {R+\K) n + P n (R-\K) r 



AR 2 -K 2 



\_Pn& mri(Pr 2^ r ) Pm^-nri^Rr 2"^" r )] 



PR 

1 2 



[e mri (R + \K) r {R + \K) n - e nri (R - \K) r {R - \K) m ] z£ 



4R 2 - K 2 



P • .R 

^mn(4-R -K" ) ~ 2~ (^n^mrs "I - i^m^-nrs)RrK t 



8m2m B 



PR 

2m 2 m^ 
4P 2 - K - 

32m2m^ 



PR 

Am2m 2 B 
) 4P 2 - 

QAm2m 2 B ' 



(B9) 



^(3,0)1 

J "ran I match 



1 



2m™ 2m 2 m 2 F 



P .i'V.P^ 

c mm L t—tL-L 



1 



a 3 !^, I £ mnr .P r z'£ • P 



2rri2m E 



2 ""-03 



a 5 (^4 + ^)£ m mZ'Si(P m + P n 



1 1 

+ 



4m B , 8m 2 m|; 8mBm E m^e 7 ™ 10- 



+ a p(u> 4 + u> 4 + «4) - a s % 



+ 



1 1 

+ 



1 



q 2 ze 



4m B , Am\ 8m2m 2 E 8mBm E mBe mia 

a\r 2 s -c 2 B )C 



+ a |(w 4 + w 4 + »y + a ^w B: 



a 3 g(w 4 + W 4 + W B )e m nrRr(i'EmRm + i^nRn) 



4m B , 8m 2 m 2 E 8mBm 2 E m E e 



mia 



+ a 3 |(w 4 + w 4 + 7w#) - a 3 |tOB 2 



+ 



1 1 

+ 



3 cl 2 ze 
+ 



4r?i4 4m 3 g , 8vii2vn 2 E 8rriBm 2 E mBe mia 
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+ a 3 \(w 4 + w' A + 7w' B ) 



a ^Wb 2 - 



4g2mi<2 



* 24 



[w 4 + u>4 + 7w' B )e mnr K r (iE m K m + iT, n K n 



a 3 \w Bl ) 5 mn P ■ R + a 6 \w A 5 mn P m R 



+ 



2m Bin e 
1 1 



34, 



1 



2r«4 4m2m E Am B m E m 2 e mi 
i 

2 



--a - 4 w Bl 



\PmRn PnRr, 



+ d c,W B 5 mn {^R m K r K m R r ^£ mr iiYj 



a 2 ZE 



4m B , 8m 2 m 2 E 8mBm 2 E m B e mia 



+ a 3 \{w 4 + w' A + Aw' B ) 
1 1 1 



a 3 \w B2 



(R n £ mr i -\- R m £ nr i)K r iYji 



a 2 z E 



4m B , 8nri2m E 8mBm E rriBe mia 

yK n £ mr i -\- K m £ nr j^R r iY^i 



1 



- a 3 \{w B2 +w B3 



Am 3 



o 3 (rf-c|)C 



21 



4 \ D2 ' o3J ^_g2m\a 

+ a 3 ±(w A + w' 4 + Aw' B )e mnr K r {iE m R m - iT, n R n 
1 



+ 



q 2 ze 



Am 3 8m2m E 8mBm 2 E m2e mia 

PmRm Rm&mriKri^i 



a 3 \w Bl 



(P K - P K 



a 2 U 



m 2 



2m E 



(BIO) 



The terms on the last line do not match, but we still must add to Eqs. (B6)-(B10) the contribution 
of the diagram with the three-gluon vertex [Fig. 2(d)], which is 

= -2iK~ 2 [25„ V R -J-(k'- K) fx J u - (k + K) V J^] + ia*\8^R^ 

+ \a 2 K~ 2 [k„k u {k> ~ K) U J V + Kk'^K + k)^] (Bll) 

and no M.^ contribution. Here J M is the current of Sec. IV B. The first lattice artifact cancels the 
last line of Eq. (BIO). The second lattice artifact vanishes upon contraction with the external-gluon 
polarization vectors. 



APPENDIX C: IMPROVED GAUGE ACTION 

In this Appendix we outline how to improve the gauge action, when axis-interchange symmetry 
is given up. The improvement program is the same as for anisotropic lattices, which has been 
worked out [24] and summarized [23]. Since it has not been published, we give the main details 
here. 

Table VI lists the interactions in the Symanzik LE£, with and without axis -interchange symme- 
try. Without axis -interchange symmetry there are eight operators. Other operators can be written 
as linear combinations of the operators in the table and total derivatives. For example, previous 
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TABLE VI: Dimension-6 gauge-field interactions that could appear in the LE£. 

With axis-interchange Without axis-interchange 

^[(D.F^iD.F^)} ti[(D 4 E) ■ (D A E)} 

^tr[(A^)(A^)] 
Y. 3 ^M{D ] B k ){D J B k )] 

tr[F^F up F pp ] tr[B • (E x E)} 

tr[B ■ (B x B)] 

tr[(D p F pu )(D p F pu )] e A tr [{D ■ E){D ■ E)} e A 

tr[(£> xB)-{Dx B)] S A 

tv[(D 4 E) ■ (D x B)} 5 E 



work [17-19] used tr[(D p F pi/ )(D p F pu )}, but we find it easier to use ti[F fJiU F up F pp ]. With the 
Bianchi identity D p F pv + D p F up + D U F PP = 0, one can show that 

±tr[(ZV^)(^-)] = ^[{D^ u ){D p F pu )\ -2tr[F pu F up F pp ] + d, (CI) 

where d denotes the omission of total derivatives that make no contribution to the action. Thus, 
only two of these three operators are needed. 

Table VI is laid out in a suggestive way: operators in the right column clearly descend from 
those in the left. It is a little harder to show that there are no more [24]. When parity and charge 
conjugation are taken into account there are 10 operators with two Ds and two Es and another 10 
where the two Es are replaced with two Bs. Of these 2x6 may be eliminated in favor of total 
derivatives and others, leaving 2 x 4 = 8 of this type. Three of these may be eliminated with the 
Bianchi identities 

D-B = 0, (C2) 
D x E = D 4 B. (C3) 

One application of the second Bianchi identity is less than obvious: 

tr[(D 4 B) ■ {D A B)\ = 2tr[B • (E x E)} - tx[{D A E) ■ (D x B)} + d. (C4) 

To find Eq. (C4) one uses Eq. (C3) for one factor of D 4 B, and then integrates by parts. In the end, 
there are 5 independent operators with two Ds and two Es or two -Bs. 

In addition, there are 6 operators with one each of D A , D, E, and B; 4 may be eliminated in 
favor of total derivatives, and another may be eliminated with a Bianchi identity, leaving 1 . Finally, 
there are the two operators tr[B ■ (E x E)} and tx[B ■ (B x B)]. Thus, the total is 8, and the list 
in Table VI is complete. 

There are three redundant interactions, corresponding to the transformations in Eqs. (2.22)- 
(2.24) that only involve gauge fields. They change the LE£ by 

£ Sym Csym + a 2 ^{e A tr[(D • E){D ■ E)] + (e A + 6 A )ti[(D x B) ■ (D x B)] 

- (2e A + 5 A + 5 E ) tv[(D 4 E) ■ (D x B)} + (e A + 5 E ) tv[(D 4 E) ■ (D 4 E)]} . (C5) 

By appropriate choice of the parameters e A , S A , and 5 E , one can remove tr[(£) ■ E)(D ■ E)] 
and two of the other three induced interactions from the LE£. Below we shall see that it is most 
convenient to choose the redundant directions as shown in the last three lines of Table VI. 
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TABLE VII: Unoriented loops on the lattice, up to length 6. 



Set i Type of loop 





' 1 r\r\ i*q 1 r~\ 1 an i i oil o c 
IClllUUIal UldAJUCllCa 


(Is 


O IJclLltll _/l tlLI LI L. L IL.i 


1 i 

LI 


Kecidngies wiin temporal long sicie 


it' 


Rectangles with temporal short side 


Is 


Spatial rectangles 


21 


"Parallelograms" with two temporal sides 


2s 


Spatial "parallelograms" 


31 


Bent rectangles with temporal bend edge 


3/' 


Bent rectangles with temporal sides, but spatial bend edge 


3s 


Spatial bent rectangles 



To construct an improved gauge action, it is enough to consider the eight classes of six-link 
loops shown in Fig. I, as well as plaquettes. Generalizing from Ref. [19], we label sets of unori- 
ented loops as in Table VII. Then let 

Si=5^2Retr[l-Z7(C)] > (C6) 

where U(C) is the product of link matrices around the curve C. The gauge action is 

S D 2 F 2 = — aSi, (C7) 

where the q are chosen so that Sc,2 F 2 > and so that classical continuum limit is correct. 

The classical continuum limit is needed not only to determine the normalization of the q, but 
also to deduce which terms in the lattice action correspond to the redundant operators of the LE£. 
The classical continuum limit of the Si is easy to find with the procedure given in Ref. [19]. For 
the plaquette terms we find 

S ot = -- / tr[E ■ E] + -y— / tx[(D A E) ■ {D A E)\ + ^ f ^tr[(AA)(A.A)], (C8) 
ci s J x L2a s J x iz j x . 

3 



Sos = -- / tr[B ■ B] + I Vtr[(/J,.5 fc )(^)L (C9) 

3r K 

where a t and a s are temporal and spatial lattice spacings, respectively. Here 

J = a t a 3 s J2 = J d * x - ( C1 °) 

It is convenient to express the six-link loops through S ot and S 0s , plus further terms of order a 2 . 
The rectangles yield 

S lt = AS ot + ^ / tr[(D 4 E) ■ (D 4 E)}, (Cll) 
a s Jx 
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S lt , = 4S 0t + a t a s [ J]tr[(A^)(A^)], (C12) 

S u = 8S 0s + ^ / Vtr[(i^B fe )(^B fe )]; (C13) 

the "parallelograms" 

S 2t = 8Sot + 4^ - Aa t a s f tr[B ■ {E x E)\ - 2a t a s f tx[{D A E) ■ (D x B)] 

J X J X 

+ a t a s [ tr[(D ■ E)(D ■ E)} - a t a s [ J]tr[(A^)(A^)], (C14) 

J X J X ■ 

An 3 f n 3 f 

S 2s = 4Sos - ^ / tr[B • (B x B)] + — tr[(B x B) ■ (D x B)] 

[ J2^[(D J B k )(D J B k )}; (C15) 

Jx • /. 



a! 



Of 

and the bent rectangles 



S'w = 8S ot + a 4 a s / tr[(D • B)(B • B)] - a t a s / ^ tr[(AA)(AA)], (C16) 

«/ £ J X ^ 

S& = 8S ot + 8S 0s - 2a t a s / ^tr[(D 4 ^ ■ (B x B)], (C17) 

S 3s = 8S 0s + ^ / tr[(B x B) • (£> x B)] — — f Vtr (B^ fc )]- (C18) 

We see immediately that the bent rectangles are the only place that the redundant interactions 
appear, so one may set c 3i , c 3t >, and c 3s at will, without sacrificing on-shell improvement. Indeed, 
the bent rectangles may be completely omitted from the improved action. 

To normalize the lattice gauge action to the classical continuum limit, one must choose 

c ot + 4(c u + c lt >) + 8c 2t + 8(c 3 t + est?) = fo, (C19) 
c 0s + 8c ls + A{c 2t + c 2s ) + 8(c 3s + c 3t ,) = £ "\ (C20) 

where £o is the bare anisotropy. At the tree level £o = CL s /a t . The essence of Eqs. (C19) and (C20) 
is to trade c ot and c 0s for the bare coupling g% and the bare anisotropy £ - 

To derive on-shell improvement conditions (at the tree level), one must allow for the transfor- 
mations in Eqs. (2.23) and (2.24). We find on-shell improvement, at the tree level, when 

(C21) 
(C22) 
(C23) 
(C24) 
(C25) 
(C26) 



£o 1<: ot — 


| - 12a* - Ax s - 4(1 + ^ 2 )x t 




|-4x t -4(4 + ^)z a> 


to^u = 


— 12 ^ X ti 


£o lc u> = 




£oCi s = 




c 2t = 


c 2s = 0, 
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£o 1 c 3 t = x t ', (C27) 
tfcw = l(x s + ^ 2 xt), (C28) 
^oC 3s = x s , (C29) 

where x t , x t >, and x s are free parameters. 

In the main text of the paper, we consider isotropic lattices, but allow for the possibility that 
heavy-quark vacuum polarization requires some asymmetry in the couplings, starting at the one- 
loop level. Thus, we consider £ = 1 an d x t = x t > = x s = x and recover [19] 

cot = c 0s = § - 24x, (C30) 

cu = cw = c ls = -± + x, (C31) 

c 2 t = c 2s = 0, (C32) 

c 3 t = c zt i = c 3s = x. (C33) 

Positivity of the action requires x < 5/72 and is guaranteed if \x\ < 1/16 [19]. Beyond the tree 
level asymmetry in these couplings may indeed arise. But the full freedom of the three redundant 
directions remains, so one may still choose c 3t = x t = 0, c 3t > = x t > = 0, and c 3s = x s = 0. 



[1] C. T. H. Davies et al [HPQCD, MILC, and Fermilab Lattice Collaborations], Phys. Rev. Lett. 92, 
022001 (2004) [arXiv:hep-lat/0304004]; C. Aubin et al. [HPQCD, MILC, and UKQCD Collabora- 
tions], Phys. Rev. D 70, 031504 (2004) [arXiv:hep-lat/0405022]; C. Aubin et al. [MILC Collabora- 
tion], Phys. Rev. D 70, 114501 (2004) [arXiv:hep-lat/0407028]. 

[2] C. Aubin et al. [Fermilab Lattice, MILC, and HPQCD Collaborations], Phys. Rev. Lett. 94, 011601 
(2005) [arXiv:hep-ph/0408306]; M. Okamoto et al, Nucl. Phys. B Proc. Suppl. 140, 461 (2005) 
[arXiv:hep-lat/0409116]. 

[3] C. Aubin et al. [Fermilab Lattice, MILC, and HPQCD Collaborations], Phys. Rev. Lett. 95, 122002 

(2005) [arXiv:hep-lat/0506030]. 
[4] I. F. Allison et al. [HPQCD and Fermilab Lattice Collaborations], Phys. Rev. Lett. 94, 172001 (2005) 

[arXiv:hep-lat/0411027]. 

[5] E. Eichten, Nucl. Phys. B Proc. Suppl. 4, 170 (1988); E. Eichten and B. R. Hill, Phys. Lett. B 234, 
511 (1990). 

[6] G. P. Lepage and B. A. Thacker, Nucl. Phys. B Proc. Suppl. 4, 199 (1988); B. A. Thacker and G. P. Lep- 
age, Phys. Rev. D 43, 196 (1991). 

[7] A. X. El-Khadra, A. S. Kronfeld, and P. B. Mackenzie, Phys. Rev. D 55, 3933 (1997) [arXiv:hep- 
lat/9604004]. 

[8] K. G. Wilson, in New Phenomena in Subnuclear Physics, edited by A. Zichichi (Plenum, New York, 
1977). 

[9] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B 259, 572 (1985). 
[10] K. Symanzik, in Recent Developments in Gauge Theories, edited by G. 't Hooft et al. (Plenum, New 

York, 1980); in Mathematical Problems in Theoretical Physics, edited by R. Schrader et al. (Springer, 

New York, 1982); Nucl. Phys. B 226, 187, 205 (1983). 
[11] A. S. Kronfeld, Phys. Rev. D 62, 014505 (2000) [arXiv:hep-lat/0002008]. 

[12] J. Harada et al, Phys. Rev. D 65, 094513 (2002) [arXiv:hep-lat/0 112044]; 71, 019903(E) (2005); 65, 
094514 (2002) [arXiv:hep-lat/01 12045]. 



42 



[13] M. B. Oktay et al, Nucl. Phys. B Proc. Suppl. 119, 464 (2003) [arXiv:hep-lat/0209150]; 129, 
349 (2004) [arXiv:hep4at/0310016]; A. S. Kronfeld and M. B. Oktay, PoS LAT2006, 159 (2006) 
[arXiv:hep4at/06 10069]. 

[14] S. Aoki, Y. Kuramashi, and S. i. Tominaga, Prog. Theor. Phys. 109, 383 (2003) [arXiv:hep- 
lat/0 107009]. 

[15] N. H. Christ, M. Li, and H. W. Lin, Phys. Rev. D 76, 074505 (2007) [arXiv:hep-lat/0608006]. 

[16] A. S. Kronfeld, in At the Frontiers of Particle Physics: Handbook of QCD, Vol. 4, edited by M. 

Shifman (World Scientific, Singapore, 2002) [arXiv:hep4at/0205021]. 
[17] P. Weisz, Nucl. Phys. B 212, 1 (1983). 

[18] P. Weisz and R. Wohlert, Nucl. Phys. B 236, 397 (1984); 247, 544(E) (1984). 

[19] M. Luscher and P. Weisz, Commun. Math. Phys. 97, 59 (1985); 98, 433(E) (1985). 

[20] M. Luscher, in Fields, Strings, and Critical Phenomena, edited by E. Brezin and J. Zinn- Justin (Else- 
vier, Amsterdam, 1990); in Probing the Standard Model of Particle Interactions, edited by R. Gupta, 
A. Morel, E. DeRafael, and E David (Elsevier, Amsterdam, 1999) [arXiv:hep4at/9802029]. 

[21] D. H. Adams and W. Lee, Phys. Rev. D 77, 045010 (2008) [arXiv:0709.0781 [hep4at]]. 

[22] G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, Phys. Rev. D 46, 4052 (1992) 
[arXiv:hep-lat/9205007]. 

[23] C. Morningstar, Nucl. Phys. B Proc. Suppl. 53, 914 (1997) [arXiv:hep4at/9608019]. 

[24] M. Alford, T. Klassen, G. P. Lepage, C. Morningstar, M. Peardon, and H. Trottier (unpublished). 

[25] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [arXiv:hep-ph/9407339]; 
55, 5853(E) (1997). 

[26] K. Osterwalder and E. Seiler, Ann. Phys. 110, 440 (1978). 

[27] M. Luscher and P. Weisz, Nucl. Phys. B 240, 349 (1984). 

[28] A. S. Kronfeld and D. M. Photiadis, Phys. Rev. D 31, 2939 (1985). 

[29] M. Di PieiTO et al. [FermiQCD Collaboration], Nucl. Phys. B Proc. Suppl. 129, 832 (2004) [arXiv:hep- 

lat/03 11027]. 
[30] M. Di Pierro, private communication. 

[31] A. S. Kronfeld, Nucl. Phys. B Proc. Suppl. 129, 46 (2004) [arXiv:hep4at/0310063]. 

[32] S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 

[33] J. Harada, S. Hashimoto, A. S. Kronfeld, and T. Onogi, Phys. Rev. D 67, 014503 (2003) [arXiv:hep- 

lat/0208004]; A. X. El-Khadra, E. Gamiz, A. S. Kronfeld and M. A. Nobes, PoS LATTICE 2007, 

242 (2007) [arXiv:071 0.1437 [hep-lat]]. 
[34] C. Bernard et al. [MILC Collaboration], Phys. Rev. D 64, 054506 (2001) [arXiv:hep4at/0104002]; 

C. Aubin et al. [MILC Collaboration], Phys. Rev. D 70, 094505 (2004) [arXiv:hep-lat/0402030]. 
[35] G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 48, 2250 (1993) [arXiv:hep4at/9209022]. 
[36] M. A. Nobes and H. D. Trottier, Nucl. Phys. B Proc. Suppl. 129, 355 (2004) [arXiv:hep4at/0309086]; 

S. Aoki, Y. Kayaba and Y. Kuramashi, Nucl. Phys. B 689, 127 (2004) [arXiv:hep4at/0401030]. 
[37] E. Follana et al. [HPQCD Collaboration], Phys. Rev. D 75, 054502 (2007) [arXiv:hep4at/0610092]. 
[38] E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu [HPQCD Collaboration], Phys. Rev. Lett. 

100, 062002 (2008) [arXiv:0706.1726 [hep4at]]. 
[39] B. A. Dobrescu and A. S. Kronfeld, Phys. Rev. Lett. 100, 241802 (2008) [arXiv:0803.0512 [hep-ph]]. 



43 



